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MODELS  FOR  PROPAGATION  OF 
NORAD  ELEMENT  SETS 


FELIX  R.  HOOTS 
RONALD  L.  ROEHRICH 


General  perturbations  element  sets  generated  by  NORAD  can 
be  used  to  predict  position  and  velocity  of  Earth-orbiting 
objects.  To  do  this  one  must  be  careful  to  use  a  prediction 
method  which  is  compatible  with  the  way  in  which  the  elements 
were  generated.  Equations  for  five  compatible  models  are 
given  here  along  with  corresponding  FORTRAN  IV  computer 
code.  With  this  information  a  user  will  be  able  to  make 
satellite  predictions  which  are  completely  compatible  with 
NORAD  predictions. 
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1.  INTRODUCTION 


NORAD  maintains  general  perturbations  element  sets  on  all 
resident  space  objects.  These  element  sets  are  periodically 
refined  so  as  to  maintain  a  reasonable  prediction  capability 
on  all  space  objects.  In  turn,  these  element  sets  are  provided 
to  users.  The  purpose  of  this  report  is  to  provide  the  user 
with  a  means  of  propagating  these  element  sets  in  time  to 
obtain  a  position  and  velocity  of  the  space  object. 

The  most  important  point  to  be  noted  is  that  not  just  any 
prediction  model  will  suffice.  The  NORAD  element  sets  are 
"mean”  values  obtained  by  removing  periodic  variations  in  a 
particular  way.  In  order  to  obtain  good  predictions,  these 
periodic  variations  must  be  reconstructed  (by  the  prediction 
model)  in  exactly  the  same  way  they  were  removed  by  NORAD. 
Hence,  inputting  NORAD  element  sets  into  a  different  model 
(even  though  the  model  may  be  more  accurate  or  even  a  numerical 
integrator)  will  result  in  degraded  predictions.  The  NORAD 
element  sets  must  be  used  with  one  of  the  models  described  in 
this  report  in  order  to  retain  maximum  prediction  accuracy. 

All  space  objects  are  classified  by  NORAD  as  near-Earth 
(period  less  than  225  minutes)  or  deep-space  (period  greater 
than  or  equal  225  minutes).  Depending  on  the  period,  the 
NORAD  element  sets  are  automatically  generated  with  the  near- 
Earth  or  deep-space  model.  The  user  can  then  calculate  the 
satellite  period  and  know  which  prediction  model  to  use. 
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2. 


THE  PROPAGATION  MODELS 


Five  mathematical  models  for  prediction  of  satellite 
position  and  velocity  are  available.  The  first  of  these,  SGP , 
was  developed  by  Hilton  §  Kuhlman  (1966)  and  is  used  for 
near-Earth  satellites.  This  model  uses  a  simplification  of 
the  work  of  Kozai  (1959)  for  its  gravitational  model  and  it 
takes  the  drag  effect  on  mean  motion  as  linear  in  time. 

This  assumption  dictates  a  quadratic  variation  of  mean  anom¬ 
aly  with  time.  The  drag  effect  on  eccentricity  is  modeled 
in  such  a  way  that  perigee  height  remains  constant. 

The  second  model,  SGP4 ,  was  developed  by  Ken  Cranford 
in  1970  (see  Lane  and  Hoots  1979)  and  is  used  for  near-Earth 
satellites.  This  model  was  obtained  by  simplification  of  the 
more  extensive  analytical  theory  of  Lane  and  Cranford  (1969) 
which  uses  the  solution  of  Brouwer  (1959)  for  its  gravitational 
model  and  a  power  density  function  for  its  atmospheric  model 
(see  Lane,  et  al  1962). 

The  next  model,  SDP4,  is  an  extension  of  SGP4  to  be  used 
for  deep-space  satellites.  The  deep-space  equations  were 
developed  by  Hujsak  (1979)  and  model  the  gravitational  effects 
of  the  moon  and  sun  as  well  as  certain  sectoral  and  tesseral 
Earth  harmonics  which  are  of  particular  importance  for  half¬ 
day  and  one -day  period  orbits. 

The  SGP8  model  (see  Hoots  1980)  is  used  for  near-Earth 
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satellites  and  is  obtained  by  simplification  of  an  extensive 
analytical  theory  of  Hoots  (to  appear)  which  uses  the  same 
gravitational  and  atmospheric  models  as  Lane  and  Cranford  did 
but  integrates  the  differential  equations  in  a  much  different 
manner. 

Finally,  the  SDP8  model  is  an  extension  of  SGP8  to  be 
used  for  deep-space  satellites.  The  deep-space  effects  are 
modeled  in  SDP8  with  the  same  equations  used  in  SDP4 . 

3.  COMPATIBILITY  WITH  NORAD  ELEMENT  SETS 

The  NORAD  element  sets  are  currently  generated  with  either 
SGP4  or  SDP4  depending  on  whether  the  satellite  is  near-Earth 
or  deep-space.  For  element  sets  sent  to  external  users,  the 
value  of  mean  motion  is  altered  slightly  and  a  pseudo-drag 
term  (n/2)  is  generated.  These  changes  allow  an  SGP  user  to 
make  compatible  predictions  in  the  following  manner.  If  the 
satellite  is  near-Earth,  then  the  pseudo-drag  term  used  in 
SGP  simulates  the  drag  effect  of  the  SGP4  model.  If  the 
satellite  is  deep-space,  then  the  pseudo-drag  term  used  in  SGP 
simulates  the  deep-space  secular  effects  of  SDP4. 

For  SGP4  and  SDP4  users,  the  mean  motion  is  first  recovered 
from  its  altered  form  and  the  drag  effect  is  obtained  from  the 
SGP4  drag  term  (B*)  with  the  pseudo-drag  term  being  ignored. 

The  value  of  the  mean  motion  can  be  used  to  determine  whether 
the  satellite  is  near-Earth  or  deep-space  (and  hence  whether 
SGP4  or  SDP4  was  used  to  generate  the  element  set) .  From 
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this  information  the  user  can  decide  whether  to  use  SGP4  or 
SDP4  for  propagation  and  hence  be  assured  of  agreement  with 
NORAD  predictions. 

The  SGP8  and  SDP8  models  have  the  same  gravitational  and 
atmospheric  models  as  SGP4  and  SDP4 ,  although  the  form  of  the 
solution  equations  is  quite  different.  Additionally,  SGP8 
and  SDP8  use  a  ballistic  coefficient  (B  term)  in  the  drag 
equations  rather  than  the  B*  drag  term.  However,  compatible 
predictions  can  be  made  with  NORAD  element  sets  by  first  cal¬ 
culating  a  B  term  from  the  SGP4  B*  drag  term. 

At  the  present  time  consideration  is  being  given  to  re¬ 
placing  SGP4  and  SDP4  by  SGP8  and  SDP8  as  the  NORAD  satellite 
models.  In  such  a  case  the  new  NORAD  element  sets  would  still 
give  compatible  predictions  for  SGP,  SGP4 ,  and  SDP4  users  and, 
for  SGP8  and  SDP8  users,  would  give  agreement  with  NORAD  predic¬ 
tions. 

4.  GENERAL  PROGRAM  DESCRIPTION 

The  five  ephemeris  packages  cited  in  section  two  have  each 
been  programmed  in  FORTRAN  IV  as  stand-alone  subroutines.  They 
each  access  the  two  function  subroutines  ACTAN  and  FM0D2P  and 
the  deep-space  equations  access  the  function  subroutine  THETAG. 
The  function  subroutine  ACTAN  is  a  two  argument  (quadrant  pre¬ 
serving)  arctangent  subroutine  which  has  been  specifically 
designed  to  return  the  angle  within  the  range  of  0  to  2 tt  .  The 
function  subroutine  FM0D2P  takes  an  angle  and  returns  the  modulo 
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by  2 7T  of  that  angle.  The  function  subroutine  THETAG  calculates 
the  epoch  time  in  days  since  1950  Jan  0  .  Q..  JJTC  *  stores  this  in 
common,  and  returns  the  right  ascension  of  Greenwich  at  epoch. 

One  additional  subroutine  DEEP  is  accessed  by  SDP4  and 
SDP8  to  obtain  the  deep'-space  perturbations  to  be  added  to  the 
main  equations  of  motion. 

The  main  program  DRIVER  reads  the  input  NORAD  2-line  element 
set  in  either  G-card  internal  format  or  T-card  transmission  for¬ 
mat  and  calls  the  appropriate  ephemeris  package  as  specified  by 
the  user.  The  DRIVER  converts  the  elements  to  the  units  of  radians 
and  minutes  before  calling  the  appropriate  subroutine.  The  ephem¬ 
eris  package  returns  position  and  velocity  in  units  of  Earth  radii 
and  minutes.  These  are  converted  by  the  DRIVER  to  kilometers  and 
seconds  for  printout. 

All  physical  constants  are  contained  in  the  constants  com¬ 
mon  Cl  and  can  be  changed  through  the  data  statements  in  the 
DRIVER.  The  one  exception  is  the  physical  constants  used  only 
in  DEEP  which  are  set  in  the  data  statements  in  DEEP. 

In  the  following  sections  the  equations  and  program  listing 
are  given  for  each  ephemeris  model.  Every  effort  has  been  made 
to  maintain  a  strict  parallel  structure  between  the  equations  and 
the  computer  code. 

5.  THE  SGP  MODEL 

The  NORAD  mean  element  sets  can  be  used  for  prediction  with 
SGP.  All  symbols  not  defined  below  are  defined  in  the  list  of 


5 


symbols  in  section  twelve.  Predictions  are  made  by  first  cal 
culating  the  constants 
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where  (t  -  tQ)  is  time  since  epoch. 


Long-period  periodics  are  included  through  the  equations 
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Solve  Kepler’s  equation  for  E  +  o’ (by  iteration  to  the 
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^1  ,  sin  u  -i 
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Short-period  perturbations  are  now  included  by 
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Then  position  and  velocity  are  given  by 


and 

r  =  r  U  +  (rv)V. 

A  FORTRAN  IV  computer  code  listing  of  the  subroutine  SGP 
is  given  below. 
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u  z  =  srwrKT^r  nine- 

VX=XMX*C0SUK-C0SN0K*SINUK 

rro — 
1-3 1 

V Y=XMY*COSUK-S INN0K*SI  NUK  " 

V  Z  =  SINIK  *C0SUK 

— TTZ— 

12 


1 


113  *  POSITION  AND  VELOCITY 


~m; - 

115 

X=RK*UX 

• 

. 

11  6~ 

Y=R<*UY 

1 1  7 

Z  *R  K  *  UZ 

rr&‘ 

x  exrr=  rd~ot*ux 

_ 

1  1  9 

Y  D0T-RD0T*UY 

TTO - 

Z  DO  r  =  RDOT  *UZ 

121 

X  DOT  =  RVDOT*VX*XDOT 

IT  2 

Y  D OT  =  RV  D  0 T*~V  Y+TD  OT 

123 

Z  D0T  =  RVD0T*VZ  +  ZDOT 

TTS - 

125 

RETURN 

T2"5  EW 
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6. 


THE  SGP4  MODEL 


The  NORAD  mean  element  sets  can  be  used  for  prediction  with 
SGP4.  All  symbols  not  defined  below  are  defined  in  the  list  of 
symbols  in  section  twelve.  The  original  mean  motion  (n0M)  and 
semimajor  axis  (ao")  are  first  recovered  from  the  input  elements 


by  the  equations 


r  I®  1  2/3 
no 


3  *2_  C3  cos  io  ' 
2  ax2  (1  -  eo2)3/2 


f 1  1  r  f  2  134  £  3^ 

ao  "  ai  (1  ■  I  51  '  61  "  81  61  ^ 


3  k2  (3  cos2  iD  -  1) 

1  ^7  u  -  O  3/2 


n"  = 


 o 


1  +  6 


a"  = 


o 


1  -  6. 


For  perigee  between  98  kilometers  and  156  kilometers,  the  value 
of  the  constant  s  used  in  SGP4  is  changed  to 

s*  =  a”  (1  -  eQ)  -  s  +  a£ 

For  perigee  below  98  kilometers,  the  value  of  s  is  changed  to 
s*  =  20/XKMPER  +  aE  . 


14 


If  the  value  of  s  is  changed,  then  the  value  of  (qQ  -s)  must 
be  replaced  by 


fq0  -  .*)' 


[(q0  -  s)4]  1/4  *  S  -  s*]4. 


Then  calculate  the  constants  (using  the  appropriate  values  of 
s  and  (qQ  -  S)4) 

0  =  cos  i 


K  - 


i"  - 


B„  ■  Cl  -  ^ 


n  =  a"e  £ 
o 


C,  *  (q„  -  s)4e4n"  (1  -  nV  7/2  taQ  (1  *  f  n2  -  4e0n 


3  3  ^2^ 

+  V  )  +  T 


(l  -  n  ) 


(  -  \  +  |  e2)  (8  +  24n2  +  3n4) ] 


C,  =  B*C, 


c„  = 


c„  = 


C%  -  s^4^A3.0n"oaE  Sin  1 


2"o  (q0  -s)V  a^602  (l-l2)'7'2 


[2n  (l  +  eQn) 


♦  i  %  +  b  3]  ‘ 


2M  2 

2  j-  [3  (1  -  39  j  (1 


a”  (1  -  n  ) 


+  1 


3  n2  .  ^  .  1  eon3}  +  |  (1  .  02)  (2n2  -  eQn 


-  eon  )  cos  2o3o] 
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C5  =  2  (q0  -  s)4C4a^602  (1  -  n2)'7/2  U  (n  *  eQ) 


*  V  1 


D2  =  4a"eCj2 


D3  -  |  aJJC2  (17aj  *  a)  Cj3 


D4  -  |  a”C5  (221a£  *  31s)  c/ 


The  secular  effects  of  atmospheric  drag  and  gravitation  are 
included  through  the  equations 

3k,(-l*  302)  3k,2  (13  -  7862  ♦  13704)  „  ,  _  , 

M-.  =  M„  ♦  [  1  *  — -  ♦  — - 3—7 - 1  o  u  V 

D  °  2ao  60  16ao  6o 


3k2  (1  -  502)  3k22  (7  -  1146 2  *  39564) 

Wrir  —  W  ^  4  +  4  ft 

Ur  O  1  Ha"  ° 


2a"  B 
o  o 


16ao  6o 


.  5k.  (3  -  360 2  *  4 9 6 4 ) 

— - -  ]"?  (t  '  V 


4a” 
o  o 


3k-2  (49  -  1903)  5k « 0  (3  -702) 

+  [ - S j  +  — - r— a -  + - 7 — g - ]  n”  (t  -t  ) 

DF  °  a"20o4  2a^  4S08  2a^4B08 


3k29 


6w  =  B*C  ,  (cos  a)  )  (t  -  t  ) 


6M 


=  (qo  -s)4B*sV!  [  (1  +  n  COS  MnP)3  -  (l  +  n  COS  Mn)3  ] 


eot1 


DF- 


16 


M  =  MnTr  +  6a)  +  6M 
p  Dr 


0)  = 


^DF  ~  ^ 


91  nMk?0  ? 

^  =  ^DF  2  °2g  2  C1  (t  "  to) 

ao  po 


e  - 


a  = 


e0  -  B*C4  (t  -  t0)  -  B*CS  (sin  Mp  -  sin  MQ) 

ao  11  -  C1  -  V  ‘  D2  (t  *  V2-  D3  (t  '  V3 


-  d4  (t  -  to)4  ]2 


IL  = 


=  M  +  a)  +  fi  +  n”  [  j  Cj  (t  -  tQ)L  +  (D2  +  2C1  )  (t  -  t0) 


♦  j  (3D3  +  120^2  +  lOCj3)  (t  -  tQ)4  +  |  (3D4  +  120^3 
+  6D22  +  30C12D2  +  15C14)  (t  -  tQ)5  ] 


=  V(1  -  e2) 


n  =  ke  /a 


/■ 


3/2 


where  (t  -  t  )  is  time  since  epoch. 


It  should  be  noted  that 
when  epoch  perigee  height  is  less  than  220  kilometers,  the  equations 


for  a  and  L  are  truncated  after  the  term,  and  the  terms  involving 
C^,  6w,  and  6M  are  dropped. 

Add  the  long-period  periodic  terms 


a  XI  =  e  cos  to 
xN 

A,  n  sin  i 

IL.  =  - 2 - ? - -  (e  cos  co)  (f  +  by) 

L  8k2a62  1  8 


3  +  50- 


17 


o 


3.0  * 


>NL  4k2afT 


ILj.  *  L  ♦ 


ayN  1  e  sin  *  *  ayNL  ' 


Solve  Kepler's  equation  for  (E  +  u>)  by  defining 


U  *=  lt  -  fl 

and  using  the  iteration  equation 

(E  +  w)i  +  1  *  (E  +  w)i  +  A (E  +  w)i 


with 

U  -  ayN  cos  (E  +  ca)i  +  axN  sin  (E  +  w)i  -  (E  +  u>) 
A(E  +  w)i  *  -  ayNsin  (E  +  a)Ji  -  axN  cos  CE  +  +  1 

and 

(E  +  w)1  =  U. 


The  following  equations  are  used  to  calculate  preliminary 
quantities  needed  for  short-period  periodics. 


e  cos 

E  -  axN  cos 

(E  + 

*  ayN 

e  sin 

E  =  axN  sin 

(E  ♦ 

“)  '  ayN 

CD 

r* 

it 

taxN2  *  ayN  5 

1/2 

Pl  = 

a  Cl  -  eL2) 

(E  +  w) 
(E  +  w) 


i  s 


r  =  a  (I  -  e  cos  E) 


1  V1T  - 

r  =  k  v  —  e  sm  E 
e  r 


rf  =  k  VeT 
0 


cos  u  *  —  [cos  (E  +  w)  -  a 


a  M  (e  sin  E) 


xN 


+ 


.  1  +Vl- 


sin  u  *  —  [sin  (E  *  w) 


axN  te  sln  E)  , 
3yN  1 


-1  ,  sm  u  . 

u  =  t  an  (  -  J 

^  cos  u  J 


L2  2 

Ar  =  ^ —  (1  ~  0  )  cos  2u 


PL 


Au  =  -  — (792  -  1)  sin  2u 

V 


3k20 

Aft  =  - j  sin  2u 

2pT  L 


3k26 

Ai  =  — =-=-  sin  i  cos  2u 
~  2  o 

2Pr 


.  k?n  7 

Ar  =  -  -  (1  -0  )  sin  2u 

PL 
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Arf  = 


k~n  9  t  ? 

[(1  -  9Z)  cos  2u  -  f  (1  -  30^)] 

L 

The  short-period  periodics  are  added  to  give  the  osculating 
quantities 

vr:;~ 

rk  =  r  [1  "  T  k2  — 2 — ~  ^302  '  +  Ar 

PL 

u^  *  u  +  Au 


=  ^  +  Aft 


*k  =  *0  *  41 


r,  =  r  +  Ar 
k 


rf^  *  r£  +  Arf  . 


Then  unit  orientation  vectors  are  calculated  by 


U  =  M  sin  uk  +  N  cos  u^. 


V  *  M  cos  u,  -  N  sin  u 
—  —  k  —  k 


where 


M  = 


M 


M 


-  sin  &k  cos  i^  1 
cos  ft^  cos  i^  > 
sin  ik 


20 


N  =  < 


N 


N. 


cos  ' 

sin  >  . 

0 


Then  position  and  velocity  are  given  by 
r  =  rkU 


and 

r  =  rkU  +  (rf)kV. 

A  FORTRAN  IV  computer  code  listing  of  the  subroutine  SGP4 
is  given  below.  These  equations  contain  all  currently  antici 
pated  changes  to  the  SCC  operational  program.  These  changes 
are  scheduled  for  implementation  in  March,  1981. 


_ ] _ * _ SGP4 _ ZL.N.G.V-.  iQ. 

2.  SUBROUTINE  S G P 4 ( I F L AG , T S I NC E ) 

3  _ COMMON  /  E  1  /  X  M  Q  ^  X^iCLDElL^aME.  'IAOvXQ  ^  CUX.NOil  N.CLL2CU - . - - - 

4  1  XNDD60/8STAR,X/Y/Z/XD0T/Y90T/ZD0T/EP0CH,DS3C 

_ 5.-. _ COMMON/  C-L/CX2  /CiC4.a.E6~A,  OOMS2X^Sv-T OXfcUtO/ - - 

6*  1  XJ  3  /  X  <  E  /  X  K  M  P  E  R  /  X  M  N  P  D  A  /  A  E 

_ i _ D  0  U  3L.E P-rt-E-C-LS-L  OiJ_£  P-O-Cii^-S-SQ - - 


3 

_ 9 _ if  hfi  a fi  .fq.  n>  fin  T.g.  i nn. 

10 


ij 
1  2 
1J 
1  4 

JL5 _ 

1  6 
J Z 

1  3 
L9 

20 

2J 

22 

2J 

24 

-15. _ 

26 

JLZ _ 

23 

30 

JSJ _ 

32 


FROM  INPUT  ELEMENTS 


JLliD  S  E  1 1  M  A  J  G  ->■  A  XX-S— 0  A  0  UP  0- 


A1=(XKE/XN0) *  *  TOTH  R  D 

_ C  Q  S  I  0  =  C  Q  S  C.X.IM.C-LJ - - - — 

THETA2=C0SI 0*C0SI3 

_ X3TidM1.=  3>»XiiEJA2r  1 .  - - - - - 

EOSQ=EO*EO 

_ 3F  I AQ. 2. ?,!■,■=■£■  CLS-3 - —  -  •  ■  - 

3ETA0=SQRT (  8  E  T  A  02 ) 

_ DEI  1=1.S»r*?»X3THMl  /  TA_1  »_Al  *flFTAft»RFTAQ2J - 

A0=A1*(1.-DELl*(.5*T0THR0+DELl*C1.*134./31.*DELl))> 

_ ftFlOa1TS»rif?*X3TMMl/fAn»An».qFTAn»flFTAn?1 -  - - 

XNODP=XNO/ ( 1 .  +DEL0) 

_ AQDP  =  A0/(1  >-.D.£LO.) - - -  - 

1 _ INITIALIZATION _ _ 

JL _ FOR  PERIGEE  LESS  THAN  220  KILOMETERS/  LdE _ LS-LM  f  LAG,  _J  S__S  El  —AM  0 _ _ 

*  THE  EQUATIONS  ARE  TRUNCATED  TO  LINEAR  VARIATION  IN  SQRT  A  AND 

*  _ QUADRATIC  VARIATION  IN  MEAN  ANOMALY. _ ALSO/  -THE  ,C3.  -LE.?.L-< — I-il  E - 

*  DELTA  OMEGA  TERM/  AND  THE  DELTA  M  TERM  ARE  DROPPED. 


3J 

34  I  S  I  M  P  =  0 

3_5 _ I  F ( ( AODP* M  .-EQ) / AE )  -LT.  (  2  2  Q  .  /  X  <  M  P  E  R  t  A  =  ).)  _L.S  LAP  =  1 _ 

36 

37  _ * _ FOR  PERIGEE  9EL0N  156  K.1/  THE  VALJES.-Q-E _ 

33  *  S  AND  QOMS2T  ARE  ALTERED 

_3 l9 _ 

40  S  4  =  S 

41  _ QOMS24  =  QOMS2T _ _ _ _ — 

42  PERIGE=(A0DP*(1.-E0)-AE)*XKMPER 

43  _ IFCPERIGE  .  G  E  .  156.)  GO  TO  10 _ 

44  S4=PER IGE-78 . 

4.5 _ I  F  (PERI  3 E  _^GJ_._93  ,J_  _G0  _LQ_2 _ 

46  S  4  =  20 . 

4  7 _ 9__  JIQMS  2  4  =  .C  ( 1_20  .  -S  4)_»A  E7XKJ3JLER 1  *  »  k _ _ 

48  S4=S4/XKMP£R+AE 

49  _ 10  PINVSQ  =  1./(AODP*AODP»3ETA02*3ETA02) _ 

50  T  S I  =  1  ./ ( A0DP-S4) 

51  _ E  TA  =  A  OP  P  »  E  0  *T  S_I _ 

52*  ETASQ=ETA*ETA 

5  3  . EEJA  =  EO*ETA _ _  _  _ 

54.  PSISQ  =  ABS  (1  .-ETAS  0  ) 

5  5* _ COEF=QOMS24  *  T  S  I  »*4 _ 

56  *  COE  FI =COE F/PS I SQ* *  3  .  5 


22 


57 

58 

59 

60 

61 

62 

63 

6  A 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 


C2  =  C0EFi*XN0DP* (AODP* ( 1 . ♦ 1 . 5*ETASCUEETA* ( 4 . ♦ETA. SO ) >  >.75* 

1  Ck2*TSI/PStS0*x3THM1#(8.O.*ETASq#(8.*FTASQ) ) ) 

C1=BSTAR#C2 

SINI0=SIN(XINCL) 

A3CVK2=-XJ3/CK2*AE**3 

C3=C0EF*TSI*A30VK2*xN0Cp*AE*SINI0/E0 

X1MTH2=1.-THETA2 

C4  =  2.*XN0DP*C0EFi*AnDP*BETA02*  <ETA* 

1  (2.*.5*ETASQ)  t-EO*  (.5*2.*ETASQ)-2.*C<2*TSI/ 

2  (AODP*PSISO) * (-3.*x3THMl* ( 1  .-2  .  *EETA«.ETASG* 

3  ( 1.5-.5*EETA) )  «-.7  5*XlMTH2*  (2.*ETASG-EETA* 

4  ( 1 .♦ETASQ) ) *COS (2,#0MEgA0) ) ) 
C5=2.*C0EFl*A0DP*BETA02* ( 1 • *2 • 75* ( ETASQ+EET A ) ♦EFTA*ETASQ) 

THETA4=ThETA2*THETA? 

TEMP1=3.*C«2*PINVSQ*XN0DP 

TEMP2=TEMP1*CK2*PINV5Q 

TEMP3=1.25*CK4#PINV5Q*PlNvSQ#XN0Dp 

XMDOT=XnODP*.5*TFMP1*BETA0*x3THMW.0625*TEmP2*BFTA0* 

1  ( 13.-78. *ThETA2+137.*THETA4) 

XlM5TH=l.-5.*THETA2 

OMGDOT=-.5*TEMP1*X1M5TH*.0625*TEMP2*(7.-114.*THFTA2> 

1  395.*THeTA4) ♦TEMPS* (3.-36. *THETA2 449. *tHETAA) 

xhdoti=-tempi*cosio 

XN0D0T=XHD0T1* (,5*TFMP2*(4.-19.*THETA2) ♦2.*TEMP3* (3.- 
1  7.*THETA2) i *COSIO 

OMGCOF=B5TAR*C3#COS (OMEGAO) 
xMCOFr-TOTHRD*COEF*RSTAR*AE/EETA 
XN0DCF=3.5*BETA02#XhD0T1*C1 
T2C0F=1.5*Cl 

XLCOF=.  125*A30vK2*STNIO*  (3.>5.*cOSIO)  /  ( l.-t-COSlO) 

AYC0F=.25*A30Vk2*SlNl0 
DELMO=(l.-*-ETA*COS(XMO)  )**3 
SlNMOrSlN (XMO) 

X7THM1=7.*THETA2-1. 

I F { I  SI MP  .EQ.  1)  GO  TO  90 
ClSQ=Cl*Cl 
D2=4.*A0DP*TSI*ClSQ 
TEMP=D2*TSI*Cl/3. 

03= ( 17.*A0Dp  +  S4) *TEmP 

D4=.5*TEMP*A0DP*TSI* (221 .*A0DP43U*54) *C1 
T3C0F=D242.*C1SQ 

T  4C0F  s . 2  5* (3.*03.C1* ( 12.*D2.10.#C1SQ) ) 

T5COF  =  .2*  <3.*DA-*>12**C1*D3«>6.*D2*D2*15.*C1SQ»  ( 

1  2.*D24ClS0n 

90  I  F|_AG  =  0 

UPDATE  FOR  SECULAR  GRAVITY  AND  ATMOSPHERIC  DRAG 

100  xMDF=xMO*xMDOT*TSlNrE 

CMGADFsoMEGA040MGDOT*TSINCE 

XNoDDF=XMODEo+XNoDoT*TsINCE 

OMEGA=0MGADF 

XMP=XMDF 

TSQ=TSINCE*TSINCE 

xNOdE=XNODDF*XNOdCF*TSO 

tempa=i.-ci*tsince 


7  X 


1 


1  1  3 


TEMPE=BSTAR»C4»TSINCE 


1  1  4 
1.151 
1  1  6 

.1 1_7__ 

IIS' 

1  1  9 


JJJL 


TEMPL=T2C0F  * t so 

_ I  F_LI  S  I  «LP_.  E  Q  0.  XO_ 

DEL0MG=0MGC0F*TSINCE 

_0  EL M  =  X M C  0 F  *  <LU  .  +.EJ.A  *CJO-S.(J<H  D.£JJ-»_»_lr HELllQJL 
TEMP=DELOMG+DELM 

XMP=XM D  F  -f  T E  W P _ 


1  20 

J  2 1 _ 

1  22 
1  2  3___ 

1  2  4 
1  25 


OMEGA=OMGAD F-TSMP 

_ I C  U  8£.=  LSil*  T.S  I N  C  E _ 

TF0UR=TSINCE*TCU3E 

_ T.EMP  A  =J  E  MPAj-  0  2_*_T  S  Q-D3»TC  U8E-D4 » T F OUR 

TEMPE  =  TEMPE  +  dSTAR*C5*  CS  I  N  (X«P>  —  S  I PJ  10) 
_ TEMPI.=TEMPL  +  T3CQF  *  T  CUBE* _ 


1  26 
1  27 
1  28 
J_2_9_ 
1  30 
.L3J_ 


1  • TF0UR*(T4C0F+TSINCE*T5C0F) 

110  A  =  AOD P»TEMPA»»2 _ 

E  =E  O-TE  MPE 

_ XL=_)^1PjL0AEG  AjLXHOD.E*  XjjgD.E*.I.£21£J _ 

8ETA=SQRT (1 ,-E*E) 
jLN=JlgLEZA**1  .5 


1  32 

1_3  3 _ 

1  34 

_1  3  5. _ 

1  36 


LONG  PERIOD  PERIODIC  $._ 


_ A_>LM= .£*£  O.SjLQ-M  E.G  AJ 


TEMP=1 . / (A*8ETA*BETA) 


1  33 

1  39 

AYNL*TEMP*AYCOF 

XLT  =  XL+XLL  .  .  _  ......  . . 

1  40 

1  4  J _ 

AYN=E*SIN(OMEGA)+AYNL 

1  4  2 
143 

★ 

SOLVE  KEPLERS  EQUATION 

1  44 

1  4*5 _ 

CAPU=FM0D2P ( XLT-XNODE) 

TEMP2=CAPJ 

1  46 
147 

DO  130  1=1^10 

SINEPW=SIN( TEMP2) 

1  48 

1  49 

C0SEPW=C0S(TEMP2) 

TEMP3  =  AXN*S  INEPW  .  .  .  .  _  _  _  _ 

1  50 

1  51 

TEMP4=AYN*C0SEPW 

TEMP5=AXN*C0SEPW 

1  52 

1  53 

T£MP6=A  YN*S  I  NE  P*l 

EPW=(CAPU-TEMP4  +  TEMP3-TEMP2  )  /  ( 1  . -TRMP5- T EMP6  )  +  T  £rlP2 

1  54 

1  55 

.  130 

IF(ABS(EPW-TEMP2)  .LE.  E  6  A )  GO  TO  140 

TEMP2=EPW 

1  56 

1  57 

★ 

SHORT  PERIOD  PRELIMINARY  QUANTITIES 

1  53 

1  59 

140 

EC0SE=TEMP5+TEMP6 

1  60 

1  61 

ES INE=TEMP3-TEMP4 

E  L  S  Q  =  A  X  N  *  A  X  N  +  A  Y  N  *  A  Y  N 

162 

1  63 

TEMP=1.-ELSQ 

PL  =  A*  TEMP 

1  6<f 

165 

R  =  A  * ( 1 .-ECOSE) 

TEMPI *1 , / R 

1  66 

1  6  7. 

RD0T=XKE*SQRT(A)*ESINE*TEMP1 

RFD0T=XKE*SQRT(PL)*TEMP1 

1  68 

TEMP2=A*TEMP1 

24 


6 9 

70 

l71. _ 

72 

L7J _ 

74 
15 
:  76 

J1 _ 

73 
79 
8C 

L8J 

:  3? 

l&a_ 

!  8  4 

185 _ 

!  36 

I  8  7 _ 

I  S3 

189 _ 

i  90 

91 _ 

!  9  2 

i  93 _ 

I  94 

I  9_5 _ 

i  96 
197 
I  93 

1  99 

200 

>01 _ 

>02 

L0.3 _ 

>04 

>05 

>06 

1 0 1 _ 

>03 

2  0  9  .  _ 
210 
>11 

>12 

>  1  3 

>  1  4 
'  15 

>  1  6 


3£  T  AL.=  SIft  T  ( I  E.MPJ - 

T  E  M  P  3  = 1 ./<1  .+36T4L) 

C0SU=TEMPZ*(  COS  EP wJ -AX.U* A  YN*-£S  I NE  * TE HP 3  )- - . — 

SINU  =  TEMP2*  (SINEP.J-AYN-AXN*ESINE*TEMP3) 

_  u=ac  la.njc  sxnu^  cas.ui - - 

SIN2J=2.*SINU*COSU 

_ C0L2lL=  2^‘XJIS  IL*_C_QJL’J-J_. - - - 

TEMP=1 . /PL 

_ TEMPI  =CK?»TEMP - 

TEMP2=TEMP1 *  T  E  M  P 


*  UPDATE  FOR  SHORT  PERIODICS 


■  RK*R*  (  1  .-1  .  5  *  TEMP  2  *8STAL  *X3THM1  )  +.  6*  TEMPI  *X  1  MTH2*COS2'J 

_ UXJ=U-*2J  »  TE  MP Z » XJ-LH-M-I-*  S  I  llZU _ 

XNODEK=XNODE+1.5*TEMP2*COSIO*SIN2U 

_ XI N  C  K  =  XI  N  C.L>  1.-5  *  I-E-M  P  2  *  C  0  S 10*  S  LN 1 0  *  C  0  S  2  U- - - - - 

RDOT<=RDOT-XM*TEMPl *  X  1 M  T  H  2  *  S I  N  2  U 

_ R  F  DlQJK-=-3.E.DQJ-JlX^J  ».I-£ MP.1.».CX  T-NIH-ZaC-Q £.2 U.+J-.J  - 

* _ _ _ QR.LEN..T4LI ON  .VECIDR.S _ 

_ SINUX  =  SIj'UU(C) _ . _ — _ 

COSUK=COS (UK) 

_ 11NJ.K=LIN  (X.INJIKJ _ 

COSIK=COS (XINCK) 

_ a.  LNN.O.  K=lS  _LNJLX_NmEKJ _ _ _ 

COSNOK-COS(XNODEK) 

_ XH  X.S- S  L1NQJL*_C-Q.1LK _ 

XMY=COSNOK*COSIK 

_ UXsXMX»S  INUK+CQSNOK* CQ.SJJ.K _ 

UY=XMY*S  INUK+SINNOK*COSU< 

_  _  U  Z  =  S  I  N  I K..*  S I  N  U  K  _  _ _ 

V  X  =  X  M  X  *  C  OSUK-C  OSNOK*SI  N  U  < 

_ VY  =  X.1Y*CQSUK-S  INN. OK  *  SIN  UK _ „ _ 

VZ=SINIK*C0SJK 


*  POSITION  AND  VELOCITY 


X=RK*UX 

Y=R  K*UY _ 

Z  =R  K  *  UZ 

XDQI=RPQTK*UX+R  LLLOXKJlVJL. 

ydot=rdotk*uy*rfdot<*vy 
Z  D 0  T  =  R  D 0  T K  *  U  Z  +  R  F  D  O-T-K*  V  Z 


RETURN  _ 
END 


o  c 


7. 


THE  SDP4  MODEL 


The  NORAD  mean  element  sets  can  be  used  for  prediction  with 
SDP4 .  All  symbols  not  defined  below  are  defined  in  the  list  of 
symbols  in  section  twelve.  The  original  mean  motion  (n^")  and 
semimajor  axis  (aQM)  are  first  recovered  from  the  input  elements 
by  the  equations 


,  e  ,2/3 

ai  ■  (  5-  1 

o 


'1  = 


3  k2  (3  cos  iQ  -  1) 
(1  -  e0Z)  3/* 


=  a. 


C1  ‘  T  6l  "  6i 


134  R  31 

tt  6i  3 


3  k2  (3  cos  iQ  -  1) 

W  Cl  -  5/2 


„»  .  n° 


1  +  6 


a1’  = 


1  -  6. 


For  perigee  between  98  kilometers  and  156  kilometers, 
the  value  of  the  constant  s  used  in  SDP4  is  changed  to 


to 


s*  =  3q  Cl  -  e„)  •  s  *  aE  • 

For  perigee  below  98  kilometers,  the  value  of  s  is  changed 
s*  =  20/XKMPER  +  aE  . 
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If  the  value  of  s  is  changed,  then  the  value  of  (qQ  -  s)  must 


be  replaced  by 


■[ 


(q0  -  S*)  =|  I(qn  -  s)  ] 


4.  1/4  , 


s  -  s 


]*. 


Then  calculate  the  constants  (using  the  appropriate  values  of 
s  and  (qQ  -  s)4) 

0  =  cos  i 


£  =  ttt 


a"  -  s 
o 


B0  -  a  -  eQ2)  VI 


n  =  a"e  £ 
o  os 


“  (Q0  -  s)4?4n^  (1  -  n2)"?//2  [a”  (1  +  j  n2  +  4eQn 


+  e  n  )  + 


3^  .  3  k2? 


(l  -  n  ) 


—  c  *  \  +  7  q2)  (8  +  24n2  +  3n4)] 


C,  =  B*C- 


=  2no  *-qo  '  s^4?4ao6o2  ^  "  r'2}~1/2  [  t2n  C1  +  eQn) 


+  7  eo  +  \  ^  - 


2k2c 


a”  (l  -  n  ) 


[3  (l  -  30 2 )  (l  +  |  n2 


-  2eon  -  jeQn3)  +  J  (1  -  e2)  (2n2  -  eQn  -  eQn3)  cos  2wo] 


M  =  [1  + 


3k?  (-1  +  302)  3k-2  (13  -  7802  +  13704) 


2  3 

2a" 
o  o 


4  7 

16ao  eo 


] n; 
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03 


=  [ 


3k2  (1  -  502)  3k22  (7  -  11462  +  39564) 

’  2ao  V  16ao\8 


5k-  (3  -  366 2  +  4994) 

+  — - 4  4 - ]  n” 

4a,,4B  0 

o  o 


n 


l 


3k26 

a"2 B  4 
o  po 


n 


It 

o 


•  •  3k,2  (40  -  1903) 

A  “  «!*[— - 4— g -  + 

2ao  0o 


5k40  (3  -  70  j 

2a,,4B  8 
o  o 


1  "o 


At  this  point  SDP4  calls  the  initialization  section  of  DEEP  which 
calculates  all  initialized  quantities  needed  for  the  deep-space 
perturbations  (see  section  ten) . 

The  secular  effects  of  gravity  are  included  by 


mdf  -  M0  *  «  (t  -  t0) 

“DF  ■  “0  +  “  (t  '  to) 

nDF  ’  So  *  K  (t  ‘  V 


where  (t  -  t  )  is  time  since  epoch.  The  secular  effect  of  drag 
on  longitude  of  ascending  node  is  included  by 


ft  =  ft 


DF 


21 

T 


noV 

a"2B  ‘ 

o  o 


(t  - 


V 


Next,  SDP4  calls  the  secular  section  of  DEEP  which  adds  the 
deep-space  secular  effects  and  long-period  resonance  effects  to 
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the  six  classical  orbital  elements  Csee  section  ten). 


The  secular  effects  of  drag  are  included  in  the  remaining 
elements  by 


1 


a  '  aDS  [1  '  C1  (t  '  V1 


L  =  M 


where  aDS,  eDS,  MDg ,  wDS,  and  ftDS,  are  the  values  of  nQ ,  eQ ,  MDF ^ 


and  ft  after  deep-space  secular  and  resonance  perturbations 
Dr 

have  been  applied. 

Here  SDP4  calls  the  periodics  section  of  DEEP  which  adds 
the  deep-space  lunar  and  solar  periodics  to  the  orbital  elements 
(see  section  ten).  From  this  point  on,  it  will  be  assumed  that 
n,  e,  I,  w,  ft,  and  M  are  the  mean  motion,  eccentricity,  inclin¬ 
ation,  argument  of  perigee,  longitude  of  ascending  node,  and 
mean  anomaly  after  lunar-solar  periodics  have  been  added. 

Add  the  long-period  periodic  terms 


e  cos  a) 


V(1  -  e7) 


A7  n  sin  i 

IL  =  -Aii! - 

I  2 

.  8k2a6^ 


Of  ^  3  +  50  ^ 

-  (e  cos  w)  (  1  + - q—  ) 


sxn  i 


o 


Et  =  IL  +  Ll 


ayN  =  e  sin  a)  +  ayNL  . 


Solve  Kepler*s  equation  for  CE  +  by  defining 


U  =  lt  -  n 


and  using  the  iteration  equation 


(E  +  uj)  i  +  1  -  CE  +  w)  i  ♦  ACE  + 


with 


and 


A(E  +  oj)  - 


U  -  a^  cos  CE  *  u)i  .  +  axN  sin  C.E  +  -  (E  +  «)  . 


Vn 


sin  (E  +  w)  ^  -  axfj  cos  ( E  +  w)  ^  +  1 


(E  +  w)1  *  U. 

The  following  equations  are  used  to  calculate  preliminary 
quantities  needed  for  the  short-period  periodics. 

e  cos  E  *  axN  cos  (E  +  w]  +  ay^  sin  CE  +  w) 

e  sin  E  =  axN  sin  CE  +  w)  -  ayN  cos  CE  +  w) 

£•  2  2.  1/2 
eL  *  (axN  *  ayN  1 


PL  =  a  (i  -  eL  ) 


r  =  a  (1  -  e  cos  E) 


r  =  k  e  sin  E 

e  r 


rf  =  k 

e  — 


r 
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cos  u  =  —  [cos  (E  +  w) 


-  a 


xN 


ayN  (e  sin  E)  ^ 

1  ♦Vl  -  eT  ^ 


sin  u  »  y  [ s in  (E  +  w) 


-  a 


yN 


axN  (e  sin  E) 


+Vl“ 


-1  ,  sin  u  s 
u  =  tan  (  — —  -  ) 


COS  u 


Ar  =  (1  -  e2)  cos  2u 

ZpT 


Au 


2  C 70 2  -1)  sin  2u 


4Pi 


3k_  6 

A  ft  =  - j  sin  2u 

2p, 


3k-0 

Ai  =  — =-*■  sin  i  cos  2u 

0  2  O 

2Pl 


k-n  7 

At  =  -  Cl  ~  0  J  sin  2u 

PL 


•  k-n  7  t,  ? 

Ar£  =  -2-  [(1  -  9  )  cos  2u  r  1  [1  •  30*)] 

PL 

The  short-period  periodics  are  added  to  give  the  osculating 
quantities 


VT~ 


rk  =  r  [1  '  2  k2 


—  ( 36 2  -  1)  ]  +  Ar 


pL 
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-1 


u  +  Au 


Then  unit 


where 


uk 

=  ^  *  ^£2 

h  3  i  +  Ai 

•  •  ■  • 
rk  =  r  +  Ar 

•  •  • 
rfk  =  rf  +  Arf  . 


orientation  vectors  are  calculated  by 


U  =  M  sin  uk  +  N  cos  u^ 


V  3  M  cos  u^  -  N  sin  u^ 


M„  ■  -  sin  Q,  cos  i, 
x  k  K 


M  =s  M  =  cos  ft,  cos  i, 
y  k  k 


lMz  -  sin  ik 


'Nx  =  cos  Qj,  ^ 


N  3  < 


Ny  =  sin  Qk 


N  =  0 
z 


Position  and  velocity  are  given  by 


-  "  rk” 

and 

r  =  rkU  +  (rf)kV. 

A  FORTRAN  IV  computer  code  listing  of  the  subroutine  SDP4 
is  given  below.  These  equations  contain  all  currently  antici¬ 
pated  changes  to  the  SCC  operational  program.  These  changes 
are  scheduled  for  implementation  in  March,  1981. 
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1 

* 

SDP4 

3  NOV 

80 

- 2 - 

*3 

- SUBROUTINE  SDP4(  I  FLAG,.  TSIMCtl  ~ 

COMMON/ El  /XM0/XN0DE0,0MEGA0/E0/XINCL/XN0,XNDT20, 

5 

- 1  OdT)51)7BTTa  R,  X,Y,Z,X  DO  T,  Y  DO  T^ZDOT,  EPOCH, 

COMMON/  Cl  /CK2/CK4,E6A,.QOMS2T,$,TOTHRD, 

D  S  53 

“  ^S~ 

7 

1  XJ3/XKE/XKf1PER*XMNPDA/«E 

DOUBLE  PRECISION  EPOCH, DS50 

~s~ 

9 

IF  <1  FLAG  .EQ.  0)  GO  TO  100 

TT5  ' 

1 1 

* 

RECOVER  ORIGINAL  MEAN  MOTION  (XNODP)  AND  SEMIMAJOR 

AXIS 

( AODP) 

— rz‘ 

1  3 

* 

“FROM  INPUT  ELEMENTS 

u 

1  5 

- A  1  =  (XKE/XN0) **  TOTHRD  " 

COSIO=COS(XINCL) 

- T5 - 

1  7 

T~HETA2  =  COSIO*COSIO 

X3THM1=3.*THET A2-1 . 

T3 

1  9 

E~OS  Q  =  EO*  EO 

BET  A02  =  1 .-EOSQ 

* 

— ?n — 

21 

- 3ETA0  =  SuftT  (B E'TXCiD 

D  EL  1 =1  .5*CK2*X3THMl / (Al*A1*8ETA0*BETA02) 

~12 - 

23 

A  0  =  A1  *  (  1  .-BETTI  *C;5*T3THRD  +  DEL1 * (1 .+  1 34./81 . *D  ELI)  )T 
DEL0=1.5*CK2*X3THM1/(A0*A0*BETA0*BETA02) 

- 

25 

‘  XNODP  =  X~N 0/(1. +DE L 0 ) 

AODP=AO/ (1 .-DELO) 

- 2"5 - 

27 

★ 

INITIALIZATION 

rs - 

29 

★ 

FOR  PERIGEE  BELOW  156  KM  *  THE  VALUES  OF 

TTD - 

31 

it 

- S  AND  G  OMS  2  T  ARE  ALTERED 

- 32 - 

33 

- STT=3 

Q0MS24=Q0MS2T 

- T5 - 

35 

PER  I  GE  =  C  AOEFP  *  (  1  . - E  OT-  A  E )  *XKMPER 

IF(PERIGE  .GE.  156.)  GO  TO  10 

-  36 

37 

- S  4  =  PER  IVE-T53 

I F  (  PE  R I G  E  .GT.  98.)  GO  TO  9 

- 33 - 

39 

- S  4  =  2TT. - 

9  Q  OM  S  2  4= ( <120.-S4)*AE/XKMPER)**4 

- 4"0 - 

41 

S  4=S4/XKMPER  +  A  E 

10  PINVSQ=1./(A0DP*A0DP*BETA02*BETA02) 

- 42 - 

43 

- S  ING  =  SIN  COME'G'Ain 

COSG=COS  COMEGAO) 

- ZT5 - 

45 

- rsi  =  1./  (  A  0  D  P  -TT4") 

ETA=AODP*EO*TS  I 

- 4~5 - 

47 

E  TASJ=ETA* ETA 

EETA=EO*ETA 

49 

PSISw  =  AcTSC1  .-ETASO) 

C  OE  F  =  Q0MS24*TS  1**4 

- 53 - 

51 

C  OE  F 1=  C  0  E  F  /  P  S  I  SO  *  *  3  .  b 

C2=COEFl*XNODP*(AODP*<1.+1.5*ETASQ+EETA*(4.+ETASQ))+.75* 

- J2 - 

53 

- T — - - CKZ*TSI7Tr$TS0*x3THM1  *<8'.0.*ETASQ*  (B.  +  EfAS 0  ))) 

C  1  =BS  T  A  R  *  C  2 

- 5*4 

55 

SINIO  =  STNTXnrCL) 

A30V<2=-XJ3/CK2*AE**3 

- 53 

- X 1MTh2=17-ThETa2 

34 


4 


57 

59 

“6TT 

61 

~6T 

63 

~~&ST 

65 


C4=2.*XN0DP*C0EF1*A0DP*BETA02*(ETA* 

.  (2.+.5*ETASQ)+E0*(,5+2.*ETASQ)-2.*CK2*TSl/ 

(A0DP*PSISQ)*(-3.*X3THM1*<1.-2.*EETA+ETASQ* 
(1 .5-. 5  *  E  ETA) )  +•  .  7  5  *  X  1  M  Tri  2*  (  2  •  *  ET  AS  Q - E  E  T A  * 

(1 ,+ETASQ) ) *C0S(2.*0MEGA0) ) ) 


66 

67 

58“ 

69 


THETA4  =  THET  X2  * THETA2 
TEMP1=3.*CK2*PINVSQ*XNOOP 


t'emp2=tempi"*cT2*pinvsq 

TEMP3=1.25*CK4*PINVSQ*PINVStt*XNOOP 


nTMOOT^XNOOP*.  5  "TEMPI  *BETAO*x3thm1+.06~2  5*TEMP2*BETAO* 

1  (13  .-73 . *THET A2  +  1 37.*THETA4) 

~ XTmT  T  H  =  1  . -5T*7HETa2~ 

OMGOOT*- . 5* TEMPI *Xl M5TH+ . 062 5*TEMP2  * (7.-1 1 4.*  THET A2  + 

“1  39  5.*THETA4)  +  TEMP3*  (3.-36.* THETA2  +  49. "THETA 4T 

XH00T1=-TEMP1 *COSIO 

XNOD0T  =  XHDOT1  +  (.5*TEMP2*(4 .-1 9. *T HE T A 2 ) +2 . * T E MP3* < 3 . - ' 
1  7.*THETA2) ) "C0SI0 


XN0DCF=3.5*BETA02*XHD0T1*C1' 

T  2C0F  =  1 . 5  * C 1 

XlC0F=.125*A30VK2*SINI0*(3.+5 

AYC0F=.25*A30VK2*SINI0 


*  COSIO) / (1 ,+COS 10) 


X7THM1=7 ,*THET A2 - 1 . 

90  IFLAG=0 

CALL  DP  I N  I  T (EOSQ,SINIO,COSIO,BETAO, AO DP, THE  TATr 


1 


SING,C0SG,BETA02,XMD0T,0MGD0T,XN0D0T,XN0DP) 


62 

83 

*  UPDATE  FOR  SECULAR  GRAVITY  AND  ATMOSPHERIC  DRAG 

8  4 

85 

100  XMDF=XMO+XMDOT*TSINCE 

*  - 

8  6 

OMGADF=OMEGa0+0MGD0Y  *  T  SINC  £ 

87 

XNODDF=XNODEO+XNODOT*T SINCE 

88 

TSU=TSlMCE*TSlNC£ 

• 

89 

XNOD£=XNODDF+XNODCF*TSQ 

90 

~TTW P  A  =  r.  -  C  1  *  T  STnT E 

91 

TEMP£=BSTAR*C4*TSINCE 

97 

TEMPL=T2C0F*TSQ 

93 

X  N  =  XN0DP 

V  A 

CALL  DPSEC(XMDF,OMGADF,XNODE,EM,XINC,XN,TSINCE) 

95 

A=(XKE/XN)**TOTHRD*TEMPA**2 

9~6  ' 

E  =  E  M-  T  E  M  P  E 

97 

XMAM=XMDF+XNODP*TEMPL 

9 "8~ 

C  A  UT  DFP  E~R  (E~7XTN  CTO  tfGA D F , XNO D E,  XMA  MT 

99 

XL=XMAM+OMGADF  +  XN0DE 

TOO" 

B  t  F  A  =  SQ  K  T ( 1  ,-E  *E ) 

101 

XN=XKE/A**1.5 

TO?” 

103 

*  LONG  PERIOD  PERIODICS 

TO  4 

105 

AXN=E*COS(OMGADF) 

1  06 

TEMP=1 . / (  A  *  9  £  T  A  *  B  £T  A ) 

107 

xll  =  temp*xlcof  *axn 

k 

T08~ 

ATN  L  =  T  b  M  P  *  A  Y“C(TF 

109 

xlt=xl+xll 

TT  0~ 

ayn=e*sin(omga  df  )  +  a  y  n  l 

* 

1  1  1 

T TT 


SOLVb  KEPLERS  EWATI0N 


113 

1 1  4 

115 

CAPU-FM0.0  2P(XLT-XN0DE.» 

TEMP2=CAPU 

116 

1 1  7 

00  130  1=1/10 

S  INEP*/  =  SIN(TEMP2> 

118 

1  1  9 

C0SEPH=C0S(TEMP2) 

TEMP3=AXN*SINEPW 

120 

121 

TEMP4=AYN*COSEPW 

TEMP5=AXN*C0SEPW 

1  22 

123 

TEMP6=AYN*SINEPW 

EPW=(CAPU-TEMP4+TEMP3~TEMP2)/(1.-TEMP5-TEMP6)+TEMP2 

ITS- 

125 

1  30 

I F  (AdS(EPW-TEMP2)  .LE.  E6A>  GO  TO  l40 

T  EMP2*EPW 

1  Zb 

127 

♦ 

SHORT  PERIOD  PRELIMINARY  QUANTITIES 

T2T 

129 

140 

EC0SE=TEMP5+TEMP6 

ITU 

131 

E SINE=TEMP3-T EMP4 

ElSQ=AXN*AXN+AYN*AYN 

132 

133 

TEMP=1 ,-ELSQ 

PL=A*TEMP 

~ TT5- 
135 

R=A*(1 .-ECOSE) 

T  EMP 1 = 1  . /R 

13  6” 

137 

TOOT  =  X KX*TQ RT(  A) *E$I  NE  ★TEMPI- 
RFD0T=X<E*SQRT(PL)*TEMP1 

“”TT3 

139 

T  EM P2  =  A  *TFMfrl 

8ETAL=SQRT (TEMP) 

tso 

1S1 

~T FM~P3  =  1  .  /  Cl  .+& ETATTX- 
C  OSU  =  T EMP2* ( COSE PW-AXN  +  AYN*E  SINE* TE MP3) 

'1ST 

143 

S  IMU  =  TEriP2*(  STRFP UFTY "N - ATTJ *ESINE*TE M”P77 

U=ACTAN(SINU/COSU) 

~T4”4 

145 

S1N2J  =  2:*S1NU*IUS~0 

C0S2U  =  2.*C0SU*C0SU“1  . 

T46 - 

147 

7TMP=1T7PT 

T  EM PI =CK2*TEMP 

r«  j 

149 

“TTMP2=TTli  FT  *TE"MP 

T5D - 

151 

★ 

— UP  OA  T  E”TOR  SHORT-" FTRIOOICS 

T TT 

153 

RK=R*  (1  .-1.5  ★TTMF?*~B~E'T  A  L  ★  X  3  T  H  M  TT  +  .6*  TEMPI  *X1MTH2*C052U 
U<=U-.25*TEMP2*X7THM1*SIN2U 

T5”4 

155 

“ XTJTTD  EK  =  XN0DE  +  1  .  5”*TE”M”F  2  *  C  0  S  1 0  *  S I  N2  U- 
XINCK=XINC+1.5*TEMP2*COSIO*SINIO*COS2U 

TT5 - 

157 

RDOT<=RDOT-XN*  TEMPI *  X  1 M  1  H  2  *  S 1n2 0 

RFD0TK=RFD0T+XN*TEMP1*(X1MTH2*C0S2U+1.5*X3THM1) 

158 

159 

★ 

ORIENTATION  VECTORS 

T5D 

161 

S INU<=S I N (UK ) 

16"2 

163 

C  OSU<  =  C  0  S  tUIO 

S  INIK  =  SIN(XINCK) 

— rrs - 

165 

TO  S  I  K  =  C  U  S  C  XTNna 

S  INN0K  =  S IN(XNODEK) 

— T66~ 

167 

C  OSNOK=CaST)TNO DEK  ) 

XMX=-SINNOK*COSIK 

— T6T5 - 

X  MY=CUSNUX*L  U ST< 
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169 

"TTO" 

171 

T77 

173 

T7T 

175 

T73- 

177 

TTT 

1  79 
T8tT 

181 

T82" 

183 

-T8T 

1  85 


UX=XMX*SINUK+COSNOK*COSU< 

- UY  =  XMY*S1N(JK  +  S  TN  NUK  +  CU  b'UT 

U  Z=S INIK  *S INUK 

- VX  =  XrtX*C.OSUK-C  Ob'NOK*  SI  N UK 

VY=XMY*COSUK-SINNOK*SINUK 
- V  Z=  S  ITJ I  K  *  C  O^UX 

■* - POSITION  AND  VELOCrn 


X  =R<*lTX 

Y =R<*UY 
T^RfC*U7 

XDOT=RDOTK*UX+RFDOTK*VX 
Y  D 0  T  =  RD'O  I  K*U  Y+"RF  DUTK  *VT 

ZDOT=RDOTK*UZ+RFDOTK*VZ 


RETURN 


8.  THE  SGP8  MODEL 

The  NORAD  mean  element  sets  can  be  used  for  prediction  with 
SGP8.  All  symbols  not  defined  below  are  defined  in  the  list  of 
symbols  in  section  twelve.  The  original  mean  motion  (n^)  and 
semimajor  axis  (a^)  are  first  recovered  from  the  input  elements 
by  the  equations 


al 


6 


1 


3  k2  (3  cos2  iQ  -  1) 
1  ^7  (1  -  eo2)3/2 


1  .  .2  134  j.  3-v 

ao  =  al  (1  '  I  61  '  61  ‘  TT  61  } 

_  k?  (3  cos2  i  -  1) 

{0  *  !  4— -OTT- 


The  ballistic  coefficient  (E  term)  is  then  calculated  from  the 
B*  drag  term  by 

B  =  2B*/p0 

where 

p  =  (2.461  x  10"!‘)  XKMPER  kg/mZ/Earth  radii 
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is  a  reference  value  of  atmospheric  density. 


Then  calculate  the  constants 


B2  =  Ire2 


6  =  cos  i 


M. 


■*  n.’  k7  2 

f  -- Cl  -  30Z) 
Z  aMV 


U,  =  -  T 


•r  nMk?  2 

|  — A  Cl  -  59  J 
z  a*  B4 


n"k. 


n,  -  r  3 


O' 


a"  B 


•  ,  n"k7  7  4 

M-  =  A-  — r^r  (13  r  780Z  +  1370  ) 

2  16  a"V 


(U-  = 


,  n*'k?  2  2 

A-  - r^r  (7  -  1140Z  +  3950  ) 

16  TV 


c  n"k/i  ?  4 

+  i  — (3  -  369Z  +  4904) 
4  a"408 

,  n"k2  2 

fl7  "  7  — 3~8’  9  C4  -  190Z) 

2  2  a"V 


C  n"k4  2 

* !  -A 9  (3  • 79  5 

z  a”4B 


Z  =  n"  +  +  M2 


a)  =  o)^  + 


TQ 


21  *  22 


a'^-s 

es? 


W 


2  ,  .  2 
a  -  1  +  e 


IB  Po 


Cq„  -  s) 4  n"  a"?4  cT1  i^7 


C-  -  1-  nM  a4  C 

12  O 


D,  = 


C^“2/a"B2 


2  9  4 

D2  =  12  +  36n  +  j  n 


2  5  4 

d3  =  I5n  +  J  n 


D4  ■  5n  ♦  x  I3 


Dr  - 


-2 


Bx  =  -k2  (1  -  se"1) 


Br,  = 


-k2  Cl  -  0  ) 


B3  -  4f  Sln  1 


C2  ~  D1D3B2 
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..  * 


C3  ' 

D4D5B3 

• 

no  * 

C1  (2  +  3n2  +  20en  +  5en3  +  e2 

+  34e2n2  +  ®2.B2B1  +  ^2  C0S  2a) 

+  Cj  sin  oj) 

C4  * 

D1D7B2 

t 

C5  ’ 

D5D8B3 

D6  * 

30n  +  n3 

D7  * 

c  25  3 

5n  ♦  —  n 

D8  = 

27  2  4 

>  l  +  =£■  n  +  n 

• 

eo 

•  -  CQ  (4n  +  n3  +  5e  +  15en2 

+  e2n  +  7e  2n3  +  D^Bj. 

+  cos  2u>  +  sin  ai) 

a/a 

=  ee  a  2 

C6  ■ 

1  n 

3  n" 

5/5 

• 

=  2aM5  (C682  +  ee) 

• 

n  = 

(e  +  e  i/K)  sK 

41 


\p/\p  -  -  nn  i|> 


-2 


c  /C  *  CA  +  4  c/C  -  a/a  -7  t|»/i|» 
0  0  o 


Cl/Cl  =  n^n"  +  4  c^a  +  VCo 

2  2 

Dg  »  6n  ♦  20e  +  liien  +  68e  n 


D10  =  20n  +  5n3  +  17e  +  68en2 


=*  72n  ♦  I8n' 


D12  =  30n  +  ion' 

d13  *  5  +  n2 


D14  *  C/C  -  2  V»/Tp 


D1S  =  2  (C6  +  ee  8*2) 


D1  =  D1  ^D14  +  °15) 


D2  =  nDll 


D3  =  nDl2 


D4  =  nD13 
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D5  "  D5D14 


C 2  R  B2  CP^^3  + 


C3  -  Bj  +  D(-D^) 


•7  n"k7  ? 

«  =  ~  j  —A  Cl  -  ™  ) 

z  a*’  B 


D16  ’  D9  *  4  D10  i  *  B1  C“lD2  4  W 


+  C2  cos  2 to  +  Cj  sin  w 


+  w  fCj  cos  a)  r  2  C2  sin  2 w) 


no  *  n  VC1  +  C1D16 

e*0  =  e  CQ/C0  -  CQ  {  (4  +  3n2  +  30en 
+  ^  e2  +  2le2n2)n  +  (5  +  I5n2 
+  31en  +  14en3)  e  +  B.^  [D-LD6 
+  D1  n  (30  +  ip-  n2)]  +  b2  [d1d7 
+  Dx  n  (5  +  7-y  n2)  ]  cos  2w 

4  B3  [°5D8  4  DSn"  <¥ 
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2  • 

+  4 n  )]  sin  w  +  oj  (Cg  cos  w 
•  2  C4  sin  2  id)} 

D1?  =  n/n"  -  (n/nM)2 

C/C  -  2  (5/5  -  C6)  5/5  +  2a"5  Cj  D1?62 

•  *2  •• 

-  2  C(-  ee  ♦  e  +  ee) 

n  =  (e  +  2e  C/C)  sC  +  n  C/C 
d18  =  C/C  -  cc/o2 

D19  ■  -  (J/’IO2  (l  +n  2)  -on  i^'2 

*D1  =  °1  (D14  +  D15)  +  D1  (D18  ‘  2D19) 

+  |  D17  +  2  a2e2B‘4  +  2  ee  B'2) 

n  =  n  D.l7  +  3  e  a  +  3  ee  a  ^ 

o  3  .L  / 

-  6  (o/a)2  ♦  4  D18  .  7  Dig] 

+  n  C1/C1  +  C1  {D^g  C1/C1 

+  Dg  n  +  D10  e  +  n2  (6  +  30en 
+  68e2)  +  ne  (40  +  3 On 2  +  272en) 
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+  e2  (17  +  68  n2) 

+  B1 

[-DiD2 

+  2D1d2  +  Dx  (n*D11  +  n2 

(72 

♦  54  n2))  ]  +  b2 

[D>3 

+  2B1°3 

*  ("n  d12  *  n2 

(30 

2 

+  30n  ))]  cos  2w 

*  B3 

[(i5Di4 

*  D5  'D18  -  2D19»  D4  * 

2Vs 

•  •  45 

+  D5  (n  °13  +  X 

nn2)] 

sin  a) 

+  oj  [(7  Cg  +  4  ee 

8-2) 

(Cj  cos  u) 

-  2  C2  sin  2w)  + 

2  C3 

COS  (a) 

-  4  C2  sin  2w  - 

“  Cc3 

sin  (a) 

+  4  C2  cos  2w) ]  } 

2  » 

2n„  -  n„  n 

O  0  0 

2  •  *•* 

-  n  n 

0  0  0 

••• 

no  1 

n  Cp'2) 

0 

45 


n 

nD  py 


6  =  - 

•  D  qy 

where  all  quantities  are  epoch  values. 

The  secular  effects  of  atmospheric  drag  and  gravitation 
are  included  by 

n  *  nQM  +  nD  [1  -  (1  -  y(t  -  tQ))P] 

e  *  eQ  +  eD  f1  '  C1  '  Y(t  -  tQ) ) q] 

“  *  «0  +  “i  t  Ct  -  tQ)  ♦  j  jprr  Zx]  *  «2  (t  -  tQ) 

o 

J!  -  C!o"  ♦  Bj  [  Ct  -  tQ)  *  j  Z1]  *  «2  Ct  -  tQ) 

M  *  Mo  *  no"  (t  '  V  *  Z1  *  M1  t(t  '  V  *  I  57  Zl] 

*  M2  Ct  -  V 

where 

zi  *  5V  {(t  *  V  *  ?T^TT  tCi  -  Y(t  -  t0))P+1  -1])  . 

If  drag  is  very  small  (— less  than  1.5  x  10  /min)  then  the 

o 

secular  equations  for  n,e,  and  should  be  replaced  by 
n  =  no"  +  n  (t  -  tQ) 
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e  -  eQ"  +  e  (t  -  tQ) 

^1=  7^o  ^  "  V 

where  (t  -  t  )  is  time  since  epoch  and  where 

2  *n 

•  L  0  f  -i  _  \ 

o 

Solve  Kepler's  equation  for  E  by  using  the  iteration 
equation 


E .  . =  E.  +  AE. 
l+l  l  l 


with 


and 


AEi  = 


M  +  e  sin  E^  -  E^ 
1  -  e  cos  E^ 


E1  =  M  +  e  sin  M  +  |  e2  sin  2M  . 

The  following  equations  are  used  to  calculate  preliminary 
quantities  needed  for  the  short-period  periodics. 

a  -  4) 2/3 


B  -  Cl  -  e2)1'2 


B  sin  E 

sin  f  =  1  -  e  cos  E 
cos  E  -  e 

cos  f  =  1  -  e  cos 


u 


f  +  0) 


aB 


1  +  e  cos  f 


p'»  =  n  a  e  sin  f 
B 


(rf ) "  =  2-|JL 

1  ? 

6r  *  -i-  — [(1  -  0  )  cos  2u 


aB' 


+  3  (1  -  302)]  -  j  sin  iQ  sin  u 


fir  =  -  n  (f)2[— |  (1  -  e2)  sin  2u  +  j  sin  i 

r  aBZ  2 

k  A 

51  -  0  if  ~T~I  sin  *0  cos  2u  '  ?  r%72  e  sin 

a  p 


a^2 


6(rf)  =  -  n  (— )  Sr  +  na  (-) 


„  sin  i 
•  a%  o 


61 


6ti  -  i  -l-T  [y  (1  -  79 2 )  sin  2u  •  3  (1 

1  aV 


-502)  (£  -  M  +  e  sin  £)] 


-  x  3  ^  [sin  i  cos  u  (2  +  e  cos  f) 

4  v? 


+  i 


0' 


2  sin  iQ/2  cos  iQ/2 


e  cos  w] 


6X  -  \  [\  (1  -  60  -  702)  sin  2u 

a  p 


-  3  (1  +  20  -  502)  (f  -  M  +  e  sin  f)] 


cos  u] 
o 
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COS  (jl) 


♦  1  ^0  sln  t 


e0 

i  +  e 


-  (2  +  e  cos  f)  cos  u] 

The  short-period  periodics  are  added  to  give  the 
osculating  quantities 
r  =  r"  +  <5r 


•  •  • 
r  =  +  5r 

r£  =  (rf) "  +  6  (rf) 

y4  =  sin  iQ/2  sin  u  +  cos  u  sin  iQ/2  6u 
+  j  sin  u  cos  iQ/ 2  61 
y5  =  sin  iQ/2  cos  u  -  sin  u  sin  iQ/2  6u 
+  j  cos  u  cos  iQ/2  61 
X  =  u  +  £2  +  6X  * 


Unit  orientation  vectors  are  calculated  by 

Ux  =  2y4  (y5  sin  X  -  y4  cos  X)  +  cos  X 

uy  =  -  2y4  (y5  cos  X  +  y4  sin  X)  +  sin  X 

Uz  -  2y4  cos  1/2 

Vx  =  2y5  Sin  X  ’  y4  C°S  ^  "  Sin  ^ 

vy  -  -2y5  (y5  cos  X  +  y4  sin  X) 

V  =  2y,-  cos  1/2 

Z  O 


+  cos  X 


where 


cos  1/2  =  V1_y42  "  y52 
Position  and  velocity  are  given  by 
r  =  r  U 

r  =  r  U  +  rf  V  , 

A  FORTRAN  IV  computer  code  listing  of  the  subroutine 
SGP8  is  given  below. 
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1 

2 

— 3_ 


s.GJ?.a_ 


4 

-5- 

6 


8 

SL¬ 


ID 
-1-1_ 
1  2 

1  4 
-13- 


1  6 
XZ 


1  8 
LSL 


20 

-2J- 


22 

-23- 


24 

2i 


26 

^Z_ 


5 

9 

30 

JJ 


32 

JUL 


34 

33 _ 

36 

LZ 


38 

39_ 

40 

4J„ 

42 

A3_ 

44 

45 

46 

47 

48 

49 


SUBROUTINE  S G P 3 ( I F L AG , T S I NC c ) 

_ COMM.ON7  E-l  ZXMO-/  XaiO-D-EX^Oai  Eg  A.Q-»C3-*-X.I  XCL-»XH-Q-,-X  N-DX2-Q-# - - 

1  XMDD60*6STAR/X/Y*Z/XD0T,YD0T,Z00T,EPOCH,D 

—C-OMAXN /  C-1  /-C  K2  a- C<  4  ,E  aA^l0_tS2X,  S-^-T  G-X4i,-i-0 , - 

1  XJ3/X<E*XKMPER,X>1NPDA/AS 

—.oauaLE— p  r.ec  i  s-LO  m—  ep.q  Cri^o_s-SC - - 

DATA  RHO/ .1 569661 5/ 


S  50 


IF  (I  FLAG  . EQ.  0)  GO  TO  100 


RECOVER  ORIGINAL  MEAN  MOTION  (XN OOP)  AND  SEMIMAJOR 
fpOiM  t  XRUI~£ LEXCXT  S.  . CA U C-UCA-LE- 


A  X  I  S  (AO 

\  CT 


kP) 
;*-X — 


<  B  TERM)  FROM  INPUT  B  *  DRAG  TERM 


A1  =  ( XKE/XNO)  **T0THRD 
cas-i.-sxcLS.-Cx ln-clo - 


THETA2=COSI*COSI 
_LXtUAaXs-3.-jL.LH  £-LA2=-1_. 


E0SQ  =  t0*  EO 
_a£XAI12=  1— -  E  Q  S  0 


BETA0=S0RT ( 3ETA02 ) 

ft  p  i  1  *  1  -  S  »  r  K  ?  »  T  T  HMLl  M  /  (A1.»_A_L*_H  E-LAO  *  BCXA02J - 

A0=Al*C1.-DELl*(.5*T0THRD+DELl*(1.+134./31 
DFI_0=1-5»CIY7»TTH/,JM/(A0»A0*-AFTA  0  I.A.Q2J - 


*  DEL  1  )  )  ) 


A0DP=A0/ ( 1 .-DELO) 
JCNQ.DP.  s.X.N-0 /  C  1-UaEL-O L 


3=2.*BSTAR/RH0 


INITIALIZATION 


I  S I MP  =0 

PQ=AQDP*B  £XA Q 2 


POM2=1 . / ( P0*P0) 
SINI=SIN  CXINCUL 


SING  =  S  IN (OMEGAO) 
JLQSGxEQ S  (OMEGAQ). 


TEMP=.5*XINCL 
SINI02  =  SIN(  TEMPI. 
C0SI02=C0S( TEMP) 
-THETA  4 sTH.ET.A2  *  *2- 


U  NM  5  T  H  =  1  ,-5.*THETA2 
UNM  T  H  2  =  1 --THETA2 


A3C0F=-XJ3/C<2*At**3 

PARQJ 1  =  3  .  *  CLK  2* POM 2 *  X.N.Q  CLP _ ___ 

PARDT2=PARDT1*C<2*P0'12 
_P  A.R  DT4=1.  25  *  C  K  4  *PQM  2 *POH 2  *X NO 0P__ . 
XMDT1=.5*PARDT1*BETAO*TTHMUN 
XGDT1  =  -.5*PARDT1  »'JNM5TH _ 


50 

51 
5  2 

53 

54 

55 


XHDT1 =-PARDT1 *C0SI 

__ X  LL  D 0  T s  X  N 0  D  P  +  X M D  T  1  + _ 

2  .0625*PARDT2*BETA0*(13.-73.*THETA2+137.*TnETA4) 

_  0M_GDJ_=XGDT1  +  _  _  _ _  ... 

1  .0 62 5  *  P ARDT 2  *  <  7 .  - 1  1  4  .  *  TH £  T  A2  +  395  .  *  TH  E  T  A4  )  +P A R 0  T 4  *  (  3  . -3 c  .  * 

2  TH£TA2+49.*THETA4) 
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XN0D0T=XHDT1+ 


1Z 

53 

59' - 

60 

--6J-- .■ - 

62 

—6.3 - 

64 

-6-5 - 

66 

-6-7 - 

63 

—6.9 _ 

70 

ZJ 

72 

-73 _ 

74 

-Z-5 _ 

76 

-7.7 _ 

78 

—73 _ 

80 

._dj _ 

82 

_8  3. _ 

84 


_ 1 _ f  ,  S  «  P  A  P  ft  T  ?  ♦  (  A  -  -  1  9  .  »T  H  F  T  A  ?  )  +  2  P.A  R  ,1T  U  C  L-  -  7  _  ft  T  :-i£lA  2-1X*  CO  SI 

T  S I  =  1 ./ (PO-S) 

_ E-I A  =  E0  *Sj»-CS-I - - — 

ETA2=ETA**2 

_ &S-I-M2  =  AdS-U-^7  -(  U^E-T-A  2-X-) - 

ALPHA2=1 .♦EOSQ 

_ £ ELA=E0*  E-T-A - 

COS2G=2. *C0SG**2-1 . 

_ 0  5  =  J-S  I  *-P-S-I  M.2 - - - 

01 =05/P0 

_ 02  =  12  .  +  E  T-A2-*  C3-&-U  4-,  UE-1-A-2) - 

03=ETA2* (15.+2.5*ETA2) 

_ 0-4-=  E-I.Aa.(  5  .  +  3-^7-5-Jt  £  LA-2  ) - 

81»CK2*TTHMUN 

_ a2-=-=£X2  tUAiM.T-H-2- - 

d3=A3COF  *SI NI 

_ C0=  ^5-*  RHO  tUCU'lST-L^Xt'WXcP  *  ACDP-*  T-S-I  *-*-4~*P-S-I-12  *  *->.  5  /  S  •)  R-T-C  AL  P  H  A  2  ) 

Cl =1 . 5*XNODP*ALPHA2**2* CO 

_ C  4  =01*0  3-*  32 - - — -  —  - 

C5=D5*D4*33 

_ XN.0J-S-C2U.-C - — 

1  (2.+6TA2* (3. +  34. *EOSQ) +5.*EETA*C4.+ETA2>  +3.5*  E JSR)  + 

_ 1 _ D.U02-*.aU _ C4-*-Ci)-S2^.>CU-SJAl.GJ - 

XNDTN=XNDT /XNODP 


*  IF  DRAG  IS  VERY  SMALL/  THE  ISIMP  FLAG  IS  SET  AND  THE 

_* _ EO.UAJ  LQ N S A-R  E-LRU  N-C A-Ii:.  0  T-0.  LLN-EJtR-VAR  LAJI-ON— Ul  MEAN 

*  MOTION  AND  QUADRATIC  VARIATION  IN  MEAN  ANOMALY 


§5 - 

8  6  I  F  (AQS  (XNDTN*XMNPD  A  )  .l.T.  2.16E-3)  GO  TO  50 

_d 1 _ Q-6-5.E..T  A»1L30.+22.S»£T  A2J _ — 
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9. 


THE  SDP8  MODEL 


The  NORAD  mean  element  sets  can  be  used  for  prediction 
with  SDP8 .  All  symbols  not  defined  below  are  defined  in  the 
list  of  symbols  in  section  twelve.  The  original  mean  motion 
(n  ")  and  semimajor  axis  (aQ")  are  first  recovered  from  the 
input  elements  by  the  equations 
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The  ballistic  coefficient  (B  term)  is  then  calculated  from 
the  B*  drag  term  by 
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is  a  reference  value  of  atmospheric  density. 
Then  calculate  the  constants 
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where  all  quantities  are  epoch  values. 

At  this  point  SDP8  calls  the  initialization  section  of  DEEP 
which  calculates  all  initialized  quantities  needed  for  the 
deep-space  perturbations  (see  section  ten) . 

The  secular  effect  of  gravity  is  included  in  mean  anomaly 
by 

MDF  =  M0  +  ^  (t  -  tQ) 

and  the  secular  effects  of  gravity  and  atmospheric  drag  are 
included  in  argument  of  perigee  and  longitude  of  ascending 
node  by 
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Next,  SDP8  calls  the  secular  section  of  DEEP  which 
adds  the  deep-space  secular  effects  and  long-period 
resonance  effects  to  the  six  classical  orbital  elements 
(see  section  ten) . 

The  secular  effects  of  drag  are  included  in  the 
remaining  elements  by 
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where  nDS,  eDS,  MDS  are  the  values  of  nQ,  eQ,  MDF  after 
deep- space  secular  and  resonance  perturbations  have  been 
applied. 

Here,  SDP8  calls  the  periodics  section  of  DEEP  which 
adds  the  deep-space  lunar  and  solar  periodics  to  the  orbital 
elements  (see  section  ten).  From  this  point  on,  it  will  be 
assumed  that  n,  e,  I,  u)>  and  M  are  the  mean  motion,  ec 
centricity,  inclination,  argument  of  perigee,  longitude 
of  ascending  node,  and  mean  anomaly  after  lunar-solar 
periodics  have  been  added. 

Solve  Kepler's  equation  for  E  by  using  the  iteration 
equation 
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The  following  equations  are  used  to  calculate  preliminary 
quantities  needed  for  the  short -period  periodics. 
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2  sm  iq/2  cos  iQ/2  J 


6  A  = 


i  4-t  [i  (1  ♦  60  -  702)  sin  2u 
z  a^S  z 

-  3  (1  +  20  -  502)  (£  -  M  +  e  sin  f )  ] 


+  T  “^77  sin  *0  [1~tV  C0S  “ 


k2a6 


-  (2  +  e  cos  f)  cos  u] 

The  short-period  periodics  are  added  to  give  the 
osculating  quantities 

r  =  r"  +  6r 
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r 


r"  +  5r 


rf  =  (rf)"  +  6(rf) 
y4  =  sin  1/2  sin  u  +  cos  iQ/2  <Su 

+  Y  sin  u  cos  i Q/ 2  61 

y5  =  sin  1/2  cos  u  -  sin  u  sin  i 

+  j  cos  u  cos  i Q/ 2  61 


X  =  u  +  ft  +  5X  . 

Unit  orientation  vectors  are  calculated  by 

Ux  =  2y4  Cy5  Sin  X  '  y4  COS  X)  + 
Uy  =  -  2y4  (y5  cos  X  +  y4  sin  X) 

Uz  =  2y4  cos  1/2 

vx  =  2y5  (y5  sin  X  -  y4  cos  X)  - 
Vy  =  -  2y5  (y5  cos  X  +  y4  sin  X) 
Vz  =  2y^  cos  1/2 

where 


and  velocity  are  given  by 
I  =  rU 

r  =  rU  +  rfV  , 


/2  6u 


cos  X 
+  sin  X 

sin  X 
+  cos  X 


Position 


A  FORTRAN  IV  computer  code  listing  of  the  subroutine 


SDP8  is  given  below. 
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1  *  SOPS  14  NOV  80 

- 2 - SUBROUTINE  S6P81 1  FL  AG#  T  5  I'N^FJ  '  " 

3  C  OMMO  N/ £  1 /  X  M  0  #  XNODEO/OMEGAO#  E  0#  X  I NC  L  # XNO#  XND  T  20/ 

- 5 - 1 - XNDD60*BlTTA  R,X  ,  Y  ,  Z,X  DOT, YDOT,ZDOT, EPOCH/ D SS  0 

5  C0MM0N/C1  /CK2,C<4,E6A,Q0MS2T  *S/TOTHRD, _ _ _ 

- 6 - - - - 1 - X J3,XK  E7K  K  M  PER,XMNPD4,AE 

7  DOUBLE  PRECISION  EPOCH/DS50 _ 

- 3  DATA  RH0/.1  569  661  5  / 

9  _ _ _ _ _ 

TD - IF  (I  FLAG  .~E'u7~(n  GO  TO  1 07] 

1  1  _ _ _  _ _ _ c _ 

~Y2 - * - ft  e  CO  U  eft —  ORIGINAL  mean  MOTION  (XNODP)  AND  SEMIMAJOR  I AXIS'  (  A  0  D  p  ) 

13  *  FROM  INPUT  ELEMENTS - CALCULATE  BALL  I S  T I C  COEFF IC I  ENT 

•  T5 - * - (6  TERMT'FftOM  input  b*  drag  term 

1  5  _ _ _ _ _ _ _ _ _ _ 

~ T6  A  1=  (XKE/XNO) *  *  TOTHRD 

17  COSI=COS(XINCL> 

T  3 - - THETA2=C0SI*T0TT  ' 

19  TTHMUN  =  3.»THETA2-1  . _ 

— 2D  E  OSQ=EO*  Eo  — 

21  BETA02=1 .-EOSO 

“22 - BETAO  =  S(JR  ITBDTAT0T7 

23  DELl»1 . 5*CK2*TTHMUN/ (A1 *A 1 *8ETA0*8E  TA02) 

— 2~it - A0*A1*(1 .-DELI *(.5*T0THRD  +  DEL1*U .  +  1 34./B1 .*DEL1>  >> 

25  DEL0=1.5*CK2*TTHMUN/(A0*A0*oETA0*3ETA02> 

— 7Z - AODP  =  AO  /  ( 1  .-DELOr 

27  X NODP=XNO/ ( 1 .+ DELO) 

“28 - B  =  2  .  *  8  S  T  A  RTU  HD 

29 

— JO - * - INITIALIZATION 

31 

— J2 - PU=AUDP*Bb’l  A 02 

33  P0m2=1 . / (PO*PO) 

J4 - S  I  N  T=  S I  N  C  X  I  N  C  L  ) 

35  S  ING=SIN (OMEGAO) 

— J 6 - C  OS G=C 0  S  COM ED7T(D 

37  TEMP=. 5*XINCL 

— 3D - S  INI02  =  SrNTTkW7 

39  C0SIU2=C0S(TEMP) 

~5D - - - - HTE  T  A  4  =  T  FTETA7'**2“ 

41  UNM5TH=1 ,-5.*THETA2 

— W - tJTTM  T  H  2  =  1  .  -  T  H  ETAT2 

43  A 3C0F =-X J 3/CK2 *AE **3 

— z TA - PAHDT1=3.*CK2*P0M2*XN0DP 

45  PARDT2=PARDT1*CK2*POM2 

'  4  6 - PAR  D  T  4  =  1  .25*CK4*pD  MT*F 3OM2*XNOOP 

47  XMDT1 =. 5*PARDT 1  * B E T A  0* T T H M U N 

4  J  "  5rDD“Tl  =-.  5  *  P’A’R' )TT*DNM  5  T  H  — 

49  XHDT1  =-PARDT1*C0SI 

5D  XLLD0T=XN0DP+XMDT1+ 

5  1  2  .06 25  *  PAR DT 2*BET AO* ( 1 3.-78. * THET A2 +1 37 T HET A4) 

52  0MuDT=XoDT1+ 

53  1  .  0625*PA  RDT2*  (  7. -1  1  4  .  *  THE  TA2+395.M  HETA  4  )  +PAR  D  T4*  (  3. -36.  * 

— 5~4  2  "  THETX2  +  49.*THETA4) 

55  XN0D0T  =  XHDT  1  + 

- - 1 - <:.5*PARDT2*<4  ."4“19.  *  TH  ET'A2)  +2  .*PARDT4*<3.-7.*THETA2))*C0SI 


57 

T  SI  =  1 ./ (PO-S) 

5"3 
.5  9 

rTA  =  EO*S.*T$  1 

7  TA?  =  ETA  **2 

6 U 

61 

P  SI M2=AdS(1  ./ci.-* E T  A T)  ) 

A  LP  H  A  2  =  1  .  +  EOSQ 

Al 

63 

E  ET A=EO*ETA 

C  OS2G  =  2 . *COSG*  *2-1  . 

64 

65 

D5  =  TSI*PSIM2“ 

D 1 =D5/P0 

6o 

67 

D2  =  l2.  +  ETA2*(36.+4.5*I:TA2) 

D3=ETA2* (15.+2.5*ETA2> 

66 

69 

"  D  4  =  E  7  A  *T3T+3  . 7~5  *~ETa7  ) 

81=C<2*TTHMUN 

- 7-0 

71 

82=-CK2*UNMTh2 

B3=A3COF *SINI  _ 

77 

73 

C"0=.5*B*RHO*QOMS2T  +  XNO DP*A0DP*TSI*+4*PSIM2**3. 57s QR T ( AL  P  H  A  2 ) 

C 1=1 . 5*XNODP*ALPHA2**?*CO 

74 

75 

C4=D1*D3*B2 

C  5=D5*D4*B3 

- 76 - 

77 

XNDT=C1*( 

1  (2.+ETA2*(3.+34.*E0SQ)+5.*EETA*(4.+ETA2)+8.5*E0SQ>+ 

- T8 

79 

1  01*02*81+  C4*C0S?G+C5*SING) 

XNDTN=XNDT/XNODP 

STJ 

81 

E  DOTs-ToTHR 0  *X  ND  TN*  C1“. -EO ) 

I FLAG=0 

- 37 

83 

- CALL  OP't  NlT<C<JSQ/SlNI  ,  C£)S  l,3EtAo;A'0'DP,THEtA2,SING,COSG, 

1  8ETA02/XLLDDT,0MGDT,XN0D0T/XN0DP> 

- 8~4 

85 

*  UPDATE  FOR  SECULAR  GRAVITY  AND  ATMOSPHERIC  DRAG 

•  -8-5 - ~ 

87 

100  Z 1=.5*XNDT*TSINCE+TSINCE 

- 8-8 - 

-'89 

- Z  7  =  3".  5"*  TO  T  H  R  D  *  Z 1  ZX7IWP 

xmamof=xmo+xlloot*tsimce 

TO 

91 

1TMGASM  =  OMEGAO  +  OMGD I  * (  SINCE  +  Z7*XG0T1 
XNODES=XNOOEO+XNODOT*TSINCE+Z7*XHOT1 

-  <72- 

93 

XN  =  'X7nT0DP 

CALL  DPSECCX MAMDF/OMG A SM/X NODES/EM/ XI NC/XN/TSINCE) 

- T? - 

95 

Xf«l=XN  +  XN  DT*  TS  I NC  t 

EM=EM+EDOT*TSINCE 

<n 

97 

HAH '^YPTA  M  0  F  ♦  Z 1  +  Z  7  +  X7TDTT- 

CALL  DPPER(EM/XINC/OMGASM/XNODES/XMAM) 

- T8 - 

99 

“7TM  A  PT=  F  M  0  D  2  P  C  X7IAM1 

— rtro - 

101 

T - SULV'b"  KbT»LbRS  bUUATrCTTO 

TO? - 

103 

rCT^TMATiT+EM*  SXNTX  MAMT*  (1.  +  EM*C0SIXMAMJ) 

00  130  1=1/10 

— nrc - 

105 

- STNT=  S  1  N  (  Z  C  2  ) 

COSE=COS (ZC2) 

— rcrs - 

107 

- Z  C5  =  1  ./Cl  .  -  b  M+'LTTS'El 

CAPE=(XMAM+EM*SINE-ZC2)* 

— 10~8 - 

109 

1  ZTS  +  Z  C  2 

I F(A8S(CAPE-ZC2)  .LE.  E6  A )  GO  TO  140 

— rro - 

1.1 1 

- 130  Z  C 2- CAFE - 

— m - 

“s - SHUKT  PtRIOD  PRELIMINARY  QUANTITIES  T- 
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113  _ * _ 

114' 

1 1  5 

140  AM  = (XKE/.XN ) *  *  TOT  HR D  •  -  .  v 

BETA2M*1 ,-EM*EM 

'  TT6 

1 1  7 

S  INOS  =  S IN(OMGASM) 

COSOS=COSCOMGASM)  *  * 

“  1 T8 

1 1  9 

A  XNM=EM* COSOS 

AY.MM=Ei<t*SIN0S 

IZ'J 

121 

PM=Art*bE  TA2M 

G 1=1 ./PM 

1~2"2 

1  2  3 

GT=T5  *C  <  2  *G  1 

G  3  =  G2  *  G  1 

rzx 

125 

6ETA  =  SuRTTBETA  2.^1 ) 

G4=.25*A3C0F*SINI 

1  2  "6 

127 

G^=.25*A3COF*Gl  ' 

SNF=BETA  *S INE*  ZC  5 

T  2  8 

1  29 

CSF=(C0SE-EM)*ZC5 

FM=ACTAN(SNF/CSF) 

T3”0 

131 

S NFG=SNF* COSOS +CSF*SIN0S 

CSFG=CSF*C0S0S-SNF*SIN0S 

1T2 

1  3  3 

SN2F2G=2 .*snfg*csfg 

CS2F2G  =  2.*CSFG**2-1  . 

13  4 
135 

EC0SF=EM*CSF 

G10=FM-XMAM+EM*SNF 

TT8T 

137 

RM=PH/( 1 .+EC0SF) 

A0\/R  =  AM/  RM 

138“ 

139 

- G1  3  =  XN*  A0\7R 

G 1 4=-Gl 3  * AOVR 

r4U~ 

141 

-  _"0"ft'=GT*TU’N M T RTMrST F 2 G - 3  ,*YTHMUN)-G4*S NTG 

DIWC=3.*G3*SINI*CS2F2G-G5*AYNfl 

‘  T  4~2 
14  3 

- D  I  =  0  I  ul  C  *  C  0S~I 

SINI2=SIN(.5*XINC> 

— T4TT- 

145 

*  UPDATE  FOR  SHORT  PERIOD  PERIODICS 

~ TTfiT- 

14  7 

SNI2DU=S  I  N  I  0  2  *  ( 

'  T4  8 

1  4  9 

1  - 5T*Tr5~*  (  1  . -7."*ThETA  2)  *  S  N2  F  2G-3  .  *  UNM1)  TH  *G“TCT  ) HTI  *  C  SF  G  *  (2.+ 

2  EC0SF))-.5*G5*THETA2*AXNM/C0SI02 

'  150 

1  5  1 

- x  LAMd  =  F  M +  0HG  AS  M+ Sf  NODES +6i*  (.  5*<  1  .  +6  .  *  CO  S  I  -  7  .  *  T  H  ET  A  2  )  *SN2>2G-r3.< 

1  <UNM5TH+2.*COSI)*GlO)+G5*SINI*<COSI*AXNM/<1.+COSI>-<2. 

r\i  no 

Ln  lo 

(r-  r- 

2 -  +  FC0S  F  )  *  CS'FG") 

Y4=SINI2*SNFG+CSFG*SNI2DU+.5*SNFG*C0SI02*DI 

— TY4~ 
155 

'  '  ~  Y  y=TTNT  2  *CSFG-SMFG*$'Nl2DU  +  .5*CSfG*C0SI02*Dt 

R  =RM+DR 

1  5'6— 

157 

- Rb0T  =  'xN'*AM*EM*SNF/BETA+G14*(2.*G2*UNMTH2*SNZF2G+G4*CSFG) 

RVD0T=XN*AM**2*8ET  A/RM  + 

TTS- 

159 

1  G  T  4*T> R  +  A  M*Gl3*SlNI*DlW ! 

16D 

1  6  1 

*  '  URTE^TATTO N  VECTORS 

162 

163 

- S  NU  ArtB=  SI  N  (  X  L'A'MBl 

CSLAMB=C0S(XLAM8) 

T64 

1  6  5 

-  - TFPTF=2'7<rCY5  *  S'- YV*  C  SL  AM3  )  “ 

UX=Y4*TEMP+CSL AMB 

T66~ 

167 

.  '  VX  =  Y  5  *  TEMP-  S  NL  A  M  9 

TEMP=2.*(Y5*CSLAMB+Y4*SNLAMB) 

T6TT 

- U  Y  =  -Y4*  r  EMP  +  SNLATTB  ~ 

-i 


169 

VY=-Y5*TEMP+CSLAM8 

T7n - 

3  7.1 

- rkMP  =  2.*SQRT(1  .  -  Y4 *TS  -  Y  5*  Y  3 ) 

UZ=Y4*TEMP 

• 

- 

172 

173 

VZ  =  Y5*TE"fTP 

17% 

175 

★ 

- P'0"S  IT  I  ON  "  ATTD  VELOCri'Y 

TT6 - 

177 

. 

— x  =r^j'x 

Y»R*UY 

‘  1 78 

1  79 

Z=R*l)Z 

XOOT»ROOT*UX+RVDOT*VX 

— rst) 

181 

- YDIH  =  KD0I*UY+R  \nrur*in 

ZDOT=RDOT*UZ+RVDOT*VZ 

— nr2 - 

183 

R  ETURN 

TffC - ETNI7 
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10.  THE  DEEP -SPACE  SUBROUTINE 

The  two  deep-space  models,  SDP4  and  SDP8,  both  access 
the  subroutine  DEEP  to  obtain  the  deep-space  perturbations  of 
the  six  classical  orbital  elements.  The  perturbation 
equations  are  quite  extensive  and  will  not  be  repeated  here. 
Rather,  this  section  will  concentrate  on  a  general  description 
of  the  flow  between  the  main  program  and  the  deep-space  sub¬ 
routines.  A  specific  listing  of  the  equations  is  available 
in  Hujsak  (1979)  or  Hujsak  and  Hoots  (1977). 

The  first  time  the  deep-space  subroutine  is  accessed  is 
during  the  initialization  portion  of  SDP4/SDP8  and  is  via 
the  entry  DPINIT.  Through  this  entry,  certain  constants 
already  calculated  in  SDP4/SDP8  are  passed  to  the  deep-space 
subroutine  which  in  turn  calculates  all  initialized  (time 
independent)  quantities  needed  for  prediction  in  deep  space. 
Additionally,  a  determination  is  made  and  flags  are  set  con¬ 
cerning  whether  the  orbit  is  synchronous  and  whether  the 
orbit  experiences  resonance  effects. 

The  next  access  to  the  deep-space  subroutine  occurs 
during  the  secular  update  portion  of  SDP4/SDP8  and  is  via 
the  entry  DPSEC.  Through  this  entry,  the  current  secular 
values  of  the  '’mean"  orbital  elements  are  passed  to  the 
deep-space  subroutine  which  in  turn  adds  the  appropriate 
deep-space  secular  and  long-period  resonance  effects  to  these 
mean  elements. 
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Code  Modification  for' Subroutine  t)EKP 

Following  line  6,  add  the  following  line 
COMMON /C2/DE2RA,  PI,  PI02,  TWOPI,  X3PI02 
Following  line  297,  add  the  following  A  lines: 
IF(X1NC  .GE.  0.)  GO  TO  90 
X1NC  =  -XINC 
XNODES  =  XNODES  +  PI 
OMGASM  =  OMGASM  -  PI 

Modify  line  298  to  become  statement  label  90. 


-/  • 


W  W  I 


w  r- 


7i  r.  r_t_ 32 - 


r>  F  Eg — S-5-A-C.E 


r  ^ - - - 

SU3R0UT  I  N€  DEEP  .  ft  c  0  .  X-I  .UC '  --~ro 

1  ,.1/r/?  CK  U.A  E6A  0aS2~T-/  S  ,-TO-TIAR-O-/ - 

_^OHMOM la  A  $t3,x<E,XK'lPER/X--1NPDA,A£  _ 

1  -  -  -  ‘■U.X3-5-C- 


o4T4  iet/5'  '  _ 

_ _ _ _ _  Ml  .  It1-' 

DATA 


ZCOSIS/ 
.^9.1444436-4 

Z  COSGS / 
^IQA-SSO-S/ 


L4La4C=4 
ZSINIS/ 


,415.4.534 - 

ZS  INGS  / 

_-w  9-63-3845  6/  - 
Z  S  I N  H  S  / 

0.0/ 


^ - 7nr^?5r^t^rrr2 .  ir«o  i »-»/ 

0  A ™  ° ® !  1  as.2A a&as* - 

’data  Q4  4/  G52 /I  .3014  993/1  .0508330/^ 

_D  A.TA-  G-5  4  A  4^41 03  3.9.8  A-- 
_  .  <r  «  DA/1T??/R00T32/^  • 


_ D  A.TA-  G  5  4  44-4108 143.A--  ~  E  _  6  3 . 7  3  9  3  7  9  2  E  -  7  / 

- - 

data  R00T54/  2.1  765303E-9/  _ _ _ _ - . 

_ _ nATft  t  H  ft  J./.4  *3.45  2AU-£=3-L — -  '  " 

.  _ LfAT-RANCE.— E-0  H — 0-E.E-P — S-E  AC_£  LttLAI-AL  IXA-U-O^A 

’  ,ftcn  BTPo^/Aa/aASQi^sujax-iO/ coso.-to / 

. . . . _... 

JidCR-=-IiiClAGJL££4.CiO - -  "  ’ 


EQ  =  EO 
XNQ— =— XNO-O.E- 


AQNV  =  1 . /AO 

X9ALU _ =-  X  T  MCJ - 

XMAO*XMO 

X£-L9  0-T-s  Oin  G  D  X.tJCN.0  CO  L 
SINQ  =  SIN(XMODEO) 

C.O  S  GL_=  _X.Q3.CXN  0  DCQO — 
OMEGAQ  =  OMEGAO 


initialize  lunar  solar  terms 


5  DAY=DS50+18261 .5dO 

J_F  ( f)  A  Y  .  E  Q  *£-S-£EJL) — 


jm__La_ic_ 


ESSLLll' InQMrf  «*=**“ 

( XNODCE ) 

UiimcEJ. 


STEM=DSIN 
r  TFM  =  DCO_S 


7 COSIL = . 91375 164- . 0 3  5630 9 6  *  C  T  E  M 

Z  S  IN  IL  =  S Q R T _X1.^JL£-Q-S  U..*  I  C Oi-LLi - - 

2Z  sVh  L  VTo  3  9  6  8  3  5 1  1  *  S  T  E  M  /  z  s  I  fl  I L 

ZCOSHL=SQR  .S-UAfclL*!  S.LNHL). - 

C =4 . 71  99 67  2*.  229971  50*  D  A  Y 
A  fn  s  5.  *.83  5. 15_Ul±  JLO J-9  itAi.6.3 

Z  MOL  =  FM0d2P(C-GAM) 
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The  last  access  to  the  deep-space  subroutine  occurs 
at  the  beginning  of  the  osculation  portion  (periodics  ap¬ 
plication)  of  SDP4/SDP8  and  is  via  the  entry  DPPER . 

Through  this  entry,  the  current  values  of  the  orbital 
elements  are  passed  to  the  deep-space  subroutine  which  in 
turn  adds  the  appropriate  deep-space  lunar  and  solar 
periodics  to  the  orbital  elements. 

During  initialization  the  deep-space  subroutine  calls 
the  function  subroutine  THETAG  to  obtain  the  location  of 
Greenwich  at  epoch  and  to  convert  epoch  to  minutes  since 
1950,  All  physical  constants  which  are  unique  to  the  deep- 
space  subroutine  are  set  via  data  statements  in  DEEP  rather 
than  being  passed  through  a  common, 

\ 

A  FORTRAN  IV  computer  code  listing  of  the  subroutine 
DEEP  is  given  below.  These  equations  contain  all  currentjy 
anticipated  changes  to  the  SCG  operational  program.  These 
changes  are  scheduled  for  implementation  in  March,  1981. 
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57 _ Z  X  =  .  3978541  6*STEM/ZSINII _ : _ 

53  l  Y=  ZCOSHL*CTErt  +  0.91  744 3 6 7 * Z S I N H L * S TE M 

5‘9J _ ZX  =  AC-TAMCIX/Zr  )  - - - - 

60  ZX=GAM+ZX-XNODCE 

61  _ ZCOSGL  =  CQS  .CZX) - 

6*2*  ZSINGL  =  SIN  (ZX) 

_6_3. _ Z MS  =  6 . 2  S..6Jl& 3 ZflJL+^21U.2iL‘L2ZZ.D :J.±D.A  1 - - - -  — 

64  ZM0S=FM0D2P CZMOS) 

66  *  DO  SOLAR  TERMS 

68  10  LS  =  0 

69  _ S.A3/XS_N.=J_^D2JD - 

70  ZC0SG=ZC0SGS 

71  _ ZSIN.G  =  ZXLNG_S _ 

72  ZC0SI=ZC0SIS 

_ LSXNL^ZSXN  IS - 

74  ZC0SH=C0SQ 

75  _ Z_S.im=S.L£LQ_ - 

76  CC=C1SS 

77  _ ZN?-ZJ>LS _ _ 

78  ZE=ZES 

79  _ ZMO  =  ZMO.S _ 

80  XNOI=1./XNQ 

__8J _ ASSIGN  30  TO  LS _ 1 _ 

82  20  A1=ZCOSG*ZCOSH+ZSING*ZCOSI*ZSINH 

,_8.3 _ AJ  =  -  ZSING*ZCO  S  H±Z,C_0_SX*X  C_CLS.lJLZ.S_tNiL _ 

84  A7=-ZC0SG*ZSINH+ZSlNG*ZC0SI*ZC0Sri 

8  5 _ A8  =  Z.S  I NG  *  Z-S-LN-I _ 

86  A9=ZSING*ZSINH+ZC0SG*ZC0SI*ZC0SH 

67  _ A10  =  ZC0SGXZSJ1U _ - 

83  A  2  =  C0SIQ*A7+  SINIQ*A8 

*89 _ A  4=  C0SIQ»A9t  SIN  IQ*A1  0 _ 

90  A  5  =-  SINIQ*A7+  C0SIQ*A3 

91  _ A.6*r_  S.  LNL9*.A-9  X_C  .0.  S.I_Q.  *  AJL2 _ 

92  C 

93  _ X  1  *  A 1  »CO  SOM  0  +  A  2  *  S  I  N  0*10 _ 

94  X2=A3*COSOMO+A4*S I NOMO 

9  5 _ X  3  =  *A  1  »S  INO  MO  +_A_2_*  CO  SLQflG _ _ _ 

96  X4=-A3*S IN0M0+A4*C0S0M0 

_  97 _ X  5  sAS  +  SINO  M  Q _ _ 

98  X6=A6*S I NOMO 

_99 _ X7  =  A5*C0S0M0 _ _ 

100  X8=A6*C0S0M0 

1  OJ _ C  _ _ 

102  Z31=12.*X1*X1-3.*X3*X3 

1  03 _ Z32«24.»X1»X2-6.»X3»X4 _ 

1 04  Z33=1 2.*X2*X2-3.*X4*X4 

_1  _as _ Z-I^j  .Xl81*A-1  ±AZ*  AZI±jU1-*  E-AX3 _ — 

106  Z2=6. * ( A1 * A3+A2*A4 ) +Z32 * EQSO 

1  0  7 _ Z3  =  3.*(A3»A3  +  A4»A4)+Z33»E  QSj _ 

10.8  Z11»-6.*Al*A5  +  EQSa  *(-24.*Xl*X7-6.*X3*X5) 

109  _ _ Z12*-6.»(A1»A6»A3*A5)»EQSQ  »(-24.»(X2»X7»Xl»X3)-6.»(X3»X6-»-X4»X5)  ) 

110  Z 1 3=-6. *A3* A6  +  EQS9  * ( -2 4 . *  X 2  *  X 8 -6 . *  X 4  *  X 6 ) 

-\7\ _ Z2136.»A2»A5  +  EQSd  »<24.»X1»XS-6.*X3»X7), _ -- 

112  Z22  =  6.*(A4*A5  +  A2*A6)«-EQSQ  *(24.*(X2*X5  +  X1*X6)-6.*(X4*X7  +  X3*X3>) 


72 


1  1  ^ 

7  2  3=6  r  *  A4*A6.fE_0Stl  *  (24...*  X  2  *  X  5-  6  -  *  X  4_*  XiD  . 

1  T  4 

M  5  _ 
r  1  6 

r  i  7 

2 1 =Z 1 +  Z1 +6SQ *Z31 

_..Z2  =  Z2  +  Z2  +  QSa*Z32 - - - - - 

Z3=Z3+Z3+8SQ*Z33 

S  3  =  C  C  *  X.N  0  I  ........  -  - -  • —  -  - 

I  1  8 

1  1  9 

S2=-.5*S3/RTEUSQ 

54=5 3*RTF 050  .  - 

1  1—7  —  .  - - 

1  20 
l  2  1 

S1=-1 5.*EQ*S4 

5S=Xl*X3+X2*X4  .  ... 

Lu.l  — — - 

122 

1  2  3 

S6=X2  *X3  +  X1  *  X  4 

S  7  =  x  2  *  X  4  -  X  1  *  x  3  . . . . 

1  24 
r  ?  s 

SE  =  S1 *  Z  N  *  S  5 

5T=52*Zl\l*(71  1+7131 

T  26 

T  2  7 

SL=-ZN*S3*<Z1+Z3-14.-6.*EQSQ) 

SGH=54*7N*(?  31+733-6.)  .  _  .  .  .  -  - 

T  23 

1  2  9 

SH=-ZN*S2*( Z21+Z23) 

IFCX9NCL*LT.5. 2J35.9  S  7  7£j=2J S  H  =&..Q _  .  - 

*  30 

1  3  1 

EE2=2.*S1*S6 

F  3  =  2 .*51*57  .  - 

1  32 

1  3  3 

XI2  =  2.*S2*Z1  2 

X_L3  =  2  .  *  S  2*  (  Z.1  3-Z  1  1  )  .  .... 

1  34 

1  35 

XL2=-2.*S3*Z2 

XL3--2.  *  S3*  ( Z3-Z1  )  .  -  ...  - .  - 

1  36 

1  3  7 

XL4=-2.*S3*(-21.-9.*EQSO)*ZE 
_ X.G.H.2=2  *54*  Z  32 _ 

1  3d 

1  39 _ 

XGH3=2. *S4* ( Z33-Z31 ) 

_ XGfcL4  =  -1  3^S4j*.1E  . 

1  40 

1  41 _ 

XH2=-2.*S2*Z22 

X  H  3  =  -  2  .  *  5  2  *  (  7  2  3-7  2  1  )  ....  .  _ 

1  42 

1  4  3 

GO  TO  LS 

1  44 

1  4_i _ 

*  0  0  LUNAR  TER. IS 

1  46 

1  4_7 _ 

30  SSE  =  SE 

5  5  T  =  5  I  — .  . 

1  48 

1  49 

SSL=SL 

_ 55H  =  SH/5INIQ _ _ _  ..  .  ..  . 

1  50 
t_5  1 

SSG=SGH-COS I Q*SSH 

5  F  2  =  F  F  2  . .  .....  ....  ....  .  -  .. 

152 

T  53 

S  I  2  =  X  12 

SL2  =  XL2  ..  _  .  . 

154 

1  5  5 

SGH2  =  XGH  2 

SJd 2_=  x  h2  . . . 

156 

1  5  7 

S  E  3  =E  3 

S  I  3  =  X  I  3  ..  .  .  . .  - . 

1  58 

1  5  9 _ 

S  L  3  =  X  L  3 

SGH3  =  XGH3 _ _ _ _ 

160 

J_6  J _ 

SH3=XH3 

S  L4  =  X.L.4 _  ....  . 

1  62 
t6  3 

SGH4=XGH4 

T  64 

T  65 _ 

ZCOSG=ZCOSGL 

2SING  =  ZSINGL  ... . -  -  ...  _ _ _  _ ... 

1  66 

J  67 

ZCOSI=ZCOSI L 

75TNT=75INTI  .  .........  - . - 

163 

ZC0SH=ZC03HL*C0SQ+ZSINHL*$ING 

7  ^ 


<4 


.6JL 


_L S.L MHs.UMQ*  ZXUSJil  -  C  Q 12  *  7.S  I  NiU, 


70 


ZN=ZNL 


7JL' 

72 

7  7 

C-C-*  C-1 L  -  --  - 

ZE=ZEL 

Z  M  0  -  Z  M  0  L  ....... 

74  * 

7  S 

ASSIGN  40  TO  LS 

76 

7  7 

40 

SSE  =  SSE+SE 

SjS  I  ~S  S I  +  S  I  . -  . 

73 

Z.9 

SSL=SSL+SL 

<:gg  =  sSG+SGh-C0SI9/SIwIQ*Sh  ...  .. 

80 

8  J 

SSH=SSH+SH/SINIQ 

82 

.3-3 

* 

GEOPOTENTIAL  RESONANCE  INITIALIZATION  FOR  12  HOUR  ORBITS 

34 

85 

I  RESFL=0 

I S  YNF1  =0  ------  - . - 

86 

8  7 

IFCXNQ.LT. (.0052359877). AND. XNQ.GT.t. 0034906585))  GO  TO  73 
fF  ._f  .  0».  y  Mil  t  G  T  .  ( E  -  ?  '  )  ..<  £  T  i  H  N 

88 

89 

IFIEQ.LT. 0.5)  RETURN 

IPFSei  =1 

90 

91 

EOC=EQ*EQSQ 

G?m  =-.  306-C  F0-.A41  *.440  .  .  .. 

92 

93 

IFCEU.GT. (.65)  )  GO  To  45 

G2l  1  s3.616-l  3-247*FQ  +  1  6-  290*FQSa_  ..  .  .  . 

I  94 
[  95 

G310*-19. 302  +  117. 390*EQ-223.»419*EQS9  +  156. 591  *EOC 

G  3  2  2  =  - 1  8  . 9068  +  10  9.  792  7.*  EQ-2.1.4 .63  .34  *£_aS_a.-til4-6^5_&l 6_*  E.O.C _  .  . 

96 
!  97 

G410*-41  .122+242. 694*EQ-47l  .094*EQSU+31  3.953  +  E0C 

G422=-1  46 .407+841 . 3  80  *  F  Q  -  1  6  2  9 . 01  4+FQSQ+  1033 . 4  35  *  FAC --- -  -  - 

1  9? 

1  99 

G520=-532.114+3017.977*EQ-5740*EQSQ+3703.276*E0C 

GO  TO  55 

?00 

>p~i 

45 

G21  1=-72. 099+331 .31  9*Ea-508.73S*£9SQ+266.724*E0C 

G  31  H  =  — 3 4  6 . 8  4  4 ♦ 1 5  3 2 . 85 1  *  E  Q-2 4  1  5 . 9  2  5  *  t QSQ  +  1 ?  4  6 . 1 1  3  *  F 0  C  . . 

?  0  2 

LOJL 

G322  =  -34 2.585  +  1  55  4. 908  +  EQ-2 366. 3  99* EQS3  +  1  2  1  5 .972  *EOC 
G410S-1Q52. 797  +  4758.636*^3-7193. 99  2  *  EQS  Q.+ 3  6  81  T  9  5  7  *  FO  C 

?  0  4 
>  0  5 

G4 22 *-35 81 .69  +  1 61 78.1  1  * E 9-2446 2 . 7 7* E 3S 9  +  1  24 22.5  2 *£0C 

IF  (EQ.GT.  (.71  5)  )  50  TO.  50 

>-06 

>07 

G520*1464. 74-4664. 75*EQ+3763.64*EQS 9 

GO  TO  5  5  ..  .  ... 

>08 

>09 

50  G520*-5149. 66  +  29936. 92*EQ-54 3 37. 36*EQ$Q  +  31324.56*E0C 

55  IF(EQ.GE.  (.  7)  )  jGQ  T_Q  _60  _  -  .  ...... 

>10 
>  1  1 

G 5  3 3  =  -91  9.2  2  7  7  +  49  3  3.61  *EQ-9064 .77*EQSQ  +  5542.21*£OC 

G  5  2  1  =  -822.71072  +  4568^61  73*£Q-6491  .4146*CQSQ  +  5337.S24*E  Qt  . 

>12 
>  1  3 

G  532  =  -853. 666  +  4690. 25*EQ- 3 624. 77*EQS Q +  5341  . 4*S0C 

GO  TO  65 

>  1  4 

’  1  5 _ 

60 

G533*-37995.78  +  161  6 16.52 *EQ-22 9 333.2 *EQSQ  + 10937  7.94 *EOC 

G521  =  -51752.104+21  3  91  3. 95*  FQ -  309463.1  6*E3S9+1  46  3  49 .4  2  *  £  0  C... 

>16 
>1  7 

65 

G532  =  -40023.88+1  70470. 39*E a- 24 2699 . 43*EQS9+1  1  5605. 32*E0C 
SINI2=SINIQ*S  INIQ  _  .. 

>18 

>19 

F220*.75*(1 ,+2.*COSia+COSQ2) 

F  221  *1  ^5  *  S  I  N 1  2 

>20 

>2*1* 

F 32 1=1 .87S*SINIQ+(1 . -2 . * C OS  IQ-3 . *C0SQ2) 

F  3  2  2  *  —  1 . 8  7  5  *  S  I  N I Q  *  (  1  .  +  2  .  *  C  0  S  I  a  ~  3  .  *  C  0  S.Q2_) 

>22 
m. . 

F441=35.*SInI2*F220 

F442  =  39.3750*SINI2  *SINI2  ..  ..  . . _.  __ 

>24 


F522  =  9.34375*SI:NI.)*(SINI2  +  (1.-2.*C05IQ-5.*C0SQ2) 


74 


J _ ♦  . 33^33333*  .C  ^7— J-4  -  *-£0-31.  J  »  A  .  ♦  CD  3j2_U - 

F  52  3  =  SlNIQ*<4.921  8751  2  *SIMI2*(-2.-4.*C3S-ia+10.*C0S-a2) 

*. _ -to  .-562-5-Q01  2-*  C-0  S  I  a  -  3,  *£-0-30-2  )-) - — 

F  5  4  2  =  29.531 25*SINia*(2.-S.*CO$ia+COSa2*(-1  2 . ♦ 3 . *  CO S  I  Q 

Jr _ *10  .*  C  03-12X4 - 

F543=29.53125*SINIQ*(-2.-8.*COSia+COSl2*(12.+8.*COSia-in.*COSG2)) 

._xM02.-=.ma<r.xN.a - 

a  i  n  v  2  =  a  a  n  v  *  a  a  n  v 

—I  S-MP-l  .  ~-3.-*  XJ10-2J-A-I -.‘FV-2 - 

TEMP  =  TEMPI *ROOT22 

_ 0.223-L  .=_' T.EMR*  E22Q*-G20.1 _ _ 

02211  =  TEMP*F221 *G21 1 

T  F  M  P  1-S_-IJ?-M.  B-lJsJilUM - 

TEMP  =  TEMPI *  R  00  T  3  2 

_ -03.2.1  q-s—tE-M-g-t-F  3  2 1  *  G  GLC - — - 

03222  =  TEMP*F322*G322 

_ XEMRJ-.  s_XEM P4-*AattV - 

TEMP  =  2  . ‘TEMPI *R00T44 

M4.lC_.s_JEMPJr-E-44-U-G.4-l  d. -  ■-  • 

04422  =  TEMP*F442*G422 

J-EM.P  1  .  s _ I  EM  Pi  *  A-GUiV- - - - 

TEMP  =  TEMP  1  *  R  00T  5  2 

0.5  220-=  .  J  EM.P  *  F.52  2J-G-S-2Q - - - 

05232  =  TEMP*F523*G532 

LE-M  P— =_ 2_*£EJie.l.*-R.Q  0  T  54- _ -- - 

05421  =  TEMP* F542*G521 

0-5-433 _ =  -LEMP.*  £-5  43*0-533 _ — - - 

XLAMO  =  XMA0+XN0DE0+XM00E0-THGR-THGR 

Q  F-A-CJ—S _ X  U_D.Q-L.-tX  .UQJlQ_T-t  X.U  30  (U-r-LH3T-=-T..hL0J - 

3FACT  =  BFACT  +SSL+S3H+SSH 

3Q_LQ_aC _ _ 


★ 


3XN.CHR.QUQ U_S_ _R£SDj1AN_CE  T..E.RM S  -INI  XL AI.IZAXLQ N - 


?Q _ LR  E3F  L  =  1 _ 

I  S  Y  N  F  L  =  1 

_ G  2.10  =J  .J2.+  E5LS  J  *_Cz2.^5+^iL12  5  *-£-'230 ) _ 

G3 1 0= 1 ,0+2.0*Easa 

_ G.3.  Q  Q  =  l  -.-CjLE  Q3Q  *  Xr.6.  -.QJt6.  6.32.33*.  £j23!L) _ _ _ 

F220=.75*C1.+C0Sia)*(1.+C0SIQ) 

F311=.9375*SINI  GL*i-LMLQ-*_C  XC.Q.S 1  3  ).  -  ^7_5jl21^*.C  233-2). 

F  330  =  1  . +  COS  I  Q 

_ F  330  =  1 . 3  75*F33Q* £3  3  Q  *  F  3  3  Q _ 

0EL1=3.*XNQ*XNQ*AQNV*AaNV 

_ DEL2  =  2.  *D  EJLJjl  E  230*.(L22Q_*i)32 _ 

0EL3  =  3.*  DELI *F330*G30Q*Q33*AQNV 

_ DELI  =D  EU*  F  3  11*  G31  0*  13  1  *A_aN4l _ 

FASX2=. 1 31 309C8 

_ FASX4=2. 88431 98 _ 

FASX6=. 37448087 

XLAMO=XMAO-*-XNODEOJ-OMEGAO-TH3R _ 

3FACT  =  XLL DOT+XP I OOT-THOT 

_ B  FAC  T  =  8  F  A,  C  J  ±.S  .5  L  +  SSS  +  S  3H _ 

80  XFACT=BF AC T-XMa 


C 


INITIALIZE  INTEGRATOR 


8J_ C 


*2 

83  . 

XLI=XLAMO. 

X  N I  =  X  N9  _  .  ..  ...... 

64 

85 

ATIME=0. DO 

5  T  F  P  P  =  7  ?  0 . 0  0  . . 

66  * 
67 

STEP 'J  =  -720. 00 

STEP?  =  259200.00  _  -  ...  -  -  . - 

86 

89 

RETURN 

90 

91 

* 

ENTRANCE  FOR  DEEP  SPACE  SECULAR  EFFECTS 

92 

93 

ENTRY  DPSEC(XLL,OMGASM,XNODES/E.vi,XINC,XN/.T) 

XLL  =  XLL  +  SSL»T  .  .....  .... 

94 

95 

OMGASM=OMGASM+SSG*T 

XNODES=XNODES+SSH*T  . . 

96 

J27 

EM=EO+SSE*T 

XINC  =  XINC_L_+S^L*T  .  . 

93 

99 

IFCIRESFL  .  EQ  .  0)  RETURN 

100  IF  (ATIME.EQ.O.DO)  GO  TO  17C  ..  ..  . .  .._ . .  ....  __  . 

:oo 

>01 

IF(T.GE.<0. DO). AND. ATIME.LT. (0.00))  GO  TO  170 

IF(T  .U.T.  (_Q  ,  D  0  )  .AND  «_A_T  IN  L.G  E-.XQ-._D  0  )  ) _  G.Q.  T  0  ±7  Q_ . 

>02 

>03 

105  IF(DA8S(T) .GE.DA3S (ATIME) )  GO  TO  120 

DELT  =  STEPP  _  .  _  _  ...  .. 

>04 

:05 

IF  (T.GE.O.DO)  DELT  =  STEPN 

110  ASSIGN  100  TO  IRET  _  _  _ _ _  _ 

o  o 

^sl  O 

GO  TO  160 

120  DELT=STEPN  .  . 

>08 

509 

IF  (T.GT.O.OO)  DELT  =  STEPP 

125  IF  (DABS (T-ATIME) -LT. STEPP)  GO  T  0  13.0.  .  . . . . 

510 

511 

ASSIGN  125  TO  IRET 

GO  TO  160  _  ...  ..  .  . .  .  .....  .  _. 

512 

;  1*3 

130  FT  =  T-ATIME 

ASSIGN  140  TO  IREJN _ _ _  ........ 

5 1  4 
515 

GO  TO  150 

140  XN  s  XMI+XNDOT*FTtXNOOT*FT*FT*a.5 

5 1  6 
517 

XL  =  XLI+XLDOT  *FT  ♦  X  N  D  0  T  *FT*FT*0. 5 

TEMP  =  -XNODES+THGR+T*THDT  ...  _  _ . . 

518 

5  1  9 

X  LL  =  XL-OMG ASM  +  TEMP 

IF  (ISYNFL.E9.0)  XLL  =  XL+TEMP+TEMP 

5  20 

5  2J _ 

C. 

RETURN 

522 

523 

c 

c 

DOT  TERMS  CALCULATED 

5  24 
525 

150  IF  (ISYNFL.EQ.O)  GO  TO  152 

XND0T=DEL1*SIN  (XLI-FASX2)+DEL2*SI  N  (2  .  *  (XL  I -F.A.S.X  4  ). ).  _ 

526 

52  7 

1  +DEL3*SIN  (3.*  (XLI-FASX6) ) 

XNDDT  =  DELI  *  COS  (  XJ__I  -  F_A  S_X_2J _ _ _ _ 

528 

529 

*  +2 . * DEL2*C0S (2 . * ( XLI-FASX4 ) ) 

*  +3. * DEL3*C0S ( 3. * ( XLI-FAS X6) ) 

530 

331 

GO  TO  154 

152  XOMI  =  0  M  E  G  A  3  +  0  M  G  D  T  *  A  T I M  E 

532 

3  3*4 

X20MI  =  XOMI  +  XOMI 

X2LI  =  XLI+XLI 

334 

3  35 

XNDOT  =  D2201 *SIN(X20MI+XLI-G22) 

*  +D221  1  *SIN(XLI-G22) 

5  3  6 

*  +D3210*SIN(XOMI+XLI-G32) 

76 


37. 


+  0322  2*A1U.L=JI 0 M  I  +ALLz&12L1. 


38 

39  _ 

40 
41_. 
42 
4J _ 

44 

45  _ 

4  6 
47 


*  +D4410*SIN(X20MI+x2LI“o44) 

_  _  +  D  44.2.2  +  S  LN-C  X2L  I-G 4  4  ).  .  - - 

*  +  05220*SIN(X0i1I+XLl-G52) 

_ * _ +  0  5  2  3  2  *  S-LNi  - XOM  I.+  XL  1 -Z52 )  . 


+D5421*SIN(X0MI+X2LI-354) 

+  P>  5433»SIN(-  XQiLI  +  X  2  L  I  r  3  5..4J. 


XNDDT 


D2201+C0S(X20MI+XLI-G22> 
.+  D.EE1  1  *ros  (  XI  I-G2.2J - 


43 

AJL 


* 

★ 


+D3210+C0S(X0HI+XLI-o32) 

+  0-32  2  2jlC.QS.<-  XOai_+_J( LI  -G3  27. 
+  D5220*COS(XO/ll+XLI-G52) 


_+ D.  57  3  7  ♦  C  0  S  (-XOM  I  +  XLI  -.G52J. 


50 

5J_ 


+2.*(D4410*C0S(X20MI+X2LI-G44) 
.+.0  4  4.27.+IL0S  (  X  2  L  L -G.4  47 _ 


52 
5  3. 
54 
55. 


56 

57  _ 


+D5421*C0S(X0HI+X2Ll-G54) 

*  ._+  D5.4  3  3  *  CO  S.C-XQH  I  +Jt  2  L  I  -  G.S.47  0 .... . 

154  XLDOT=XNl+X  FACT 

_ XjN DDT  ._=  _ X  N D  0 J.LX L Q.QJ _ 

TO 


GO 


IR6TN 


53 

59. 

60 

61 


..C. 

C 

C 


INTEGRATOR 


160 


ASSIGN  165 

_G.Q_JLQ _ 15..Q _ 


T  0  IRETN 


62 

63  _ 

64 

6.5 _ 

66 

.67 _ 


165  X L I  = 


XL  I +XLOOT  *  0  S  L  T  +XNDQT+STEP2 
XNI  =  X  N.L+  XNOOT  *  0  E  L  T.  +  X  N  D  OJ  *  ST.EP.2.. 
ATIHE=AT  IHE+DELT 

_ S.Q_  1 0_L R £J _ 


XPQ..CH  RES.IA.RX 


63 
6  9_ 

70 

71 

72 

LL 


1  70 


1  75 
L3.0 


74 

.7.5 _ 

7  6 
77  .  . 
73 
79 


l  F _ .( T . G  E. 0 . 0 0 ) 

DELT=STEPN 

_G  0 _ T  .0 _ L30 _ 

DELT  =  STEPP 

AT  I  HE  =  0  «_DlQ _ 

XNI =XNQ 

_xl  i=.x_lam_0 _ 

TO 


G  0_  T0  J7.5_. 


GO 


125 


__C_ 

C 

c. 


ENTRANCES  FOR  LUNAR-SOLAR  PERIOOICS 


30 

81 

32 

83 

84 

85 


ENTRY  0  P  P  E  R  (  E  Mr  X  I  N  C  /  OH  3.  A  S  M/  X  NO  0.  ES./  X  LL  )  _ 

SIN  IS  =  SIN(XINC) 


COS  IS  =  COS  (  X  I  NC)  _  _ _  .  _ 

IF  (OAQS(SAVTSN-T).LT. (30.00)) 
S  A  V  T  S  N  =  T _ 


GO  TO  210 


86 

87 

88 

89 

90 
9_1_ 
9  2 


205 


ZM=ZMOS  +  ZNS *  T 
Z  F_=  Zrt +2_._*  Z  E  S  *  S  I N  (  Z  M  )_ 
SINZF=SIN  (ZF  ) 

F  2  =  .  5  *  S.IN  Z  F  +S  INZ.F-.  25 
F  3=-. 5*SINZF*C0S  (ZF) 
_S.G5.ES  E  2  *  F  2+-S  E3-+.E3 _ 


S IS=S 12  *  F2  +  S I 3*F3 


9J 

SI SsSL?* F?+SL3*F3  +  SI 4*S LN7F  -  -  -  - 

94 

SGHS=SGH2*F2+S3H3*F3+SGH4*SINZF 

SHS=SH?»F?+SH3»F3  .  -  —  -  -  - 

96 

9  7 

ZM=Z  MOL  +  Z  NL  *  T 

7f:  7  '■T+  2  *  7  F  LftS-I  C  Z  M  )  -  —  .  - - 

93  * 

99 

S  I  N  Z  F  =  S  I  N  <  Z  F  ) 

F  ?=  -  S  #S  T  P4 7  F  *  SI  M2  F-  ?  5  _  _  .  .. 

00 
n  i 

F3  =  -.5*SINZF*COS  ( Z  F  > 

s EL  =  ££2*  E?+_£A*_E_3  -  — -  - 

02 
n  3 

SIL=XI2*F2+XI3*F3 

S  1  1  =  X  L 2  *  F  ?  +  X  1  3  *  F  3  +  X  1.  4  *  S.LM Z_f  -  . .  ..... 

04 

ns 

SGHL  =  XGH2*F2  +  XGH3*F3>XGH4*SI.NZF 

S  Ml  si  H?  *  F  ?♦  X  H  3  *  F  3  .  -  -  -  - - - - 

06 

07 

PE=SES+SEL 

PINTsSI  S  +  SI  1  -  -  -  -  — . -  - 

08 

0  9 

?  i  n 

PL=SLS+SIL 

PfSHsSfXHS  +  SGMI  -  -  -  -  • 

10 

1  1 

PH=SHS+SHL 

yTMr  a  iciMr.p.iMC -  -  ...  - 

1  2 

1  x 

EM  =  EM+PE 

TCfnnun  .  i  t  .  (  ,  2  ) )  GO  LQ  Z2C  -  - 

1  4 

1  s 

r 

GO  TO  218 

1  6 

17 

H 

! 

APPLY  PERIODICS  DIRECTLY 

1  3 

1  9 

218 

PHsPH/SINIQ 

PRHaPGH-CnS  I  3  *  P  H  .....  .  — 

20 
^2 1 

OMGASM=OMGASM+PGH 

XNnpFS  =  *Nnr>F$+PH  .  -  —  .  . 

At 

,  ?  3 

XLL  =  XLL+PL 
■SO  TO  ?30  

.24 
.  2*5 

C 

C 

APPLY  PFRTOnirS  WITH  1YDDANF  tt Q ft T F I C A T T 0 N - 

.26 

.2.7 _ 

.23 

.29 

C 

_ 2-2.0 _ S-LMflK  =  S  I N  (JLN.QlDE.S-) —  - 

C0S0K=CQS(X NODES) 

ALFDP  =  SINIS*SJ-Nd« .  . -  - 

.30 
t  31 

dETDP=SINIS*COSOK 

DAIF=PH»rOSOK+PINC»COSIS»SINQK  ..  ...  .  ... 

;32 

4  3  3 

D3ET=-PH*SIN0K+PINC*C0SIS*C0S0K 

ALFDP  =  ALFDP  +  DALF  .  .  .  .  -  — ■ 

434 

4  35 

3ETDP=BETDP+D3ET 

XLS  =  XLL+OM GAS M+C0SIS*X NODES  .  .  _.  . 

4  3  6 

4  37 

dls=pl+pgh-pinc*xnodes*sinis 

XLS  =  XLS  +  rH.S  - 

433 

4  39 

XN0DES=ACTAN(ALFDP/3ETDP) 

X  LL  =  XL L  +  PL 

4  40 

4  4  1 

230 

OMGASM  =  XLS-XLL-C0S(XINC)*XN0DES 

CONTINUE  .  ...  . - 

442 

443 

RETURN 

END  

•4 
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11.  DRIVER  AND  FUNCTION  SUBROUTINES 


The  DRIVER  controls  the  input  and  output  function  and 
the  selection  of  the  model.  The  input  consists  of  a  program 
card  which  specifies  the  model  to  be  used  and  the  output 
times  and  either  a  G-card  or  T?card  element  set. 

The  DRIVER  reads  and  converts  the  input  elements  to 
units  of  radians  and  minutes.  These  are  communicated  to 
the  prediction  model  through  the  common  El.  Values  of  the 
physical  and  mathematical  constants  are  set  and  communicated 
through  the  commons  Cl  and  C2,  respectively. 

The  program  card  indicates  the  mathematical  model  to 
be  used  and  the  start  and  stop  time  of  prediction  as  well 
as  the  increment  of  time  for  output.  These  times  are  in 
minutes  since  epoch. 

In  the  interest  of  efficiency  the  DRIVER  sets  a  flag 
(IFLAG)  the  first  time  the  model  is  called.  This  flag  tells 
the  model  to  calculate  all  initialized  (time  independent)  , 
quantities.  After  initialization,  the  model  subroutine 
turns  off  the  flag  so  that  all  subsequent  calls  only  access 
the  time  dependent  part  of  the  model.  This  mode  continues 
until  another  input  case  is  encountered, 

The  DRIVER  takes  the  output  from  the  mathematical 
model  (communicated  through  the  common  El)  and  converts  it 
to  units  of  kilometers  and  seconds  for  printout. 
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The  function  subroutine  ACTAN  is  passed  the  values 
of  sine  and  cosine  in  that  order  and  it  returns  the  angle 
in  radians  within  the  range  of  0  to  2it.  The  function  sub¬ 
routine  FM0D2P  is  passed  an  angle  in  radians  and  returns 
the  angle  in  radians  within  the  range  of  0  to  2it  .  The 
function  subroutine  THETAG  is  passed  the  epoch  time  exactly 
as  it  appears  on  the  input  element  cards.*  The  routine 
converts  this  time  to  days  since  1950  Jan  _Q ,_Q  UTC^,  stores  this 
in  the  common  El ,  and  returns  the  right  ascension  of  Greenwich  at 

epoch  (in  radians). 

FORTRAN  IV  computer  code  listings  of  the  routines 
DRIVER,  ACTAN,  FM0D2P ,  and  THETAG  are  given  below. 


*If  only  one  year  digit  is  given  (as  on  standard  G-cards) 
the  program  assumes  the  80  decade.  This  may  be  overridden 
by  putting  a  2  digit  year  in  columns  30-31  of  the  first 
G-card. 
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1 

* 

DRIVER  3  N0V 

80 

- 2 - 

3 

* 

WGS-72  PHYSICAL  AND  GEOPOTENTIAL  CONSTANTS' 

5 

* 

C<2=  .  5  *  J  2  *  A  E  *  "*”2  CK4=-.375*J4*AE**4 

- 6 

7 

DOUBLE  PRECISION  EPOCH/DS5"0 

C0MMON/EWXM0/XN0DEO/0MEGA0/E0/XINCL/XN0/XNDT20/XNDD60/3STAR/ 

4 

— s 

9 

- 1 -  X/Y/Z/XDOT/YDOT/ZDOT/EPOCH/DSbd 

COMMON/C1/CK2/CK4/E6A/QOMS2T/S/TOTHRD/ 

TO 

1  1 

1 —  X  J  J/"X"KE  /  XKMPER/XMNPDA/AE 

C0MM0N/C2/DE2RA/PI/PI02/TW0PI/X3PI02 

~T2 

1  3 

DATA  I H  G7THGT- 

DATA  DE2RA/E6A/PI/PI02/Q0/SO/TOTHRD/TW0PI/X3PI02/XJ2/XJ3/ 

1'4 

1  5 

1  "  XJ4/XKE/XKMPER/XMNPDA/AE/.1?4532925E-1/1 .£-6/ 

2  3.14159265/1 . 5 70 796 3 3 /I  20 . 0/ 78 . 0/ . 666666 67/ 

15 

1  7 

- 4 - 6718  3  f8T3z  4. 71  2  3  8893/1.08  2616  E  -  3/ -. 25388 1 E-5  / 

5  -1. 65597 E- 6/. 7436 69 161E-1/6378. 135/1440. /I./ 

T8 

1  9 

DIMENSION  I  S  E  T  ("5  ) 

CHARACTER  ABUF*80(2) 

20 

21 

- DATA  <  I  SET(  I  > /  I  =  f73T73HSGP/4HSGP4/4HSDP4,4HSGP8/4HSDP8/ 

22“ 

23 

♦ 

SELECT  EPHEMER  IS  TYPE  AND  OUTPUT  TIMES 

2  4” 

25 

C  <2=75 * x 72*  a F*  *7 

CK4=-.375*XJ4* AE**4 

26 

27 

- CTO MS7T=  (WO^ST) ~*  AT7X  KM’PETTT *  *  5* 

S  =A  E  *  (  1 ,+SO/XKMPER) 

“28 

29 

2~1TEXD  C5770tf)  I  f  PT7  TT7T F  /  D  E  LT 

I F( IEPT.LE.O)  STOP 

30 

31 

- rUETP  =  0 

% 

32  - 

33 

- * - 

- REA  D”TIT  WEAN  ELfcMbNIS  1-kOM  Z  card  I  ( TRArts)  OR  GTTFrrTR~KD  FO  fTM  A 1 

• 

« 

“T4 

35 

—  READ  (5/706T  ABUT 

D EC  ODE ( A8UF ( 1  ) /707 )  ITYPE 

36 

37 

IFCITYPE.  E"57TrTG~)  GO  TO  5 

DECODE  (AdUF/702)  E P OC H / X N D T 2 0 / XN DD 6 0 / I E X P / 8 S T A R/ I B E X P / X I N C L / 

— 33 - 

39 

- 1 - XNODEO/EO/OM'EGW/XMO/XN'O- 

GO  TO  7 

30 

41 

- 5  D  EC  OD'E  (  ABUF  /  70  1  )  E"POCH/XMO/XNODEO/OME  G  A^7E"0TxT  NCL/XN0/XNDT20/ 

1  XNDD60/IEXP/8STAR/I8EXP 

32 

43 

“7 - I  F  (  X  N  0 .  L  E  .  U  .  )  ST  OP” 

WRITE (6/704)  ABUF/ISET ( IEPT) 

- 3T5 - 

45 

- 1  F(IEPT.GT.b)  'GO'  TO  WO 

XNDD6O=XNDD6O*(10.**IEXP) 

45 

47 

X  NODEO  =  X  NOD E 0*  DE 2 R A 

0MEGA0=0MEGA0*DE2RA 

4“8 

49 

XTTO^TMO  *  '6F2 TTS- 
XINCL=XINCL*DE2RA 

- 5D - 

51 

- TEMP  =  TUO'PI/XMNPDA/XMnP'DA 

X  N0  =  XN0*  TEMP*XMNPDA 

5"2  " — 

53 

-  XND  I20=XNDT20*  1  EMP 

XNDD60=XNDD60*TEMP/XMNPDA 

54 

55 

★ 

INPUT  CHECK  FOR  PERIOD  VS  EPHEMER  IS  SELECTED 

- 515 - 

it 

- PERIOD  GE  22b  MINUTE 3  IS  DEEP  SPACE-' 

81 

57 


"5TT 

59 

*6tr 

61 


63 


■yzr 

65 

“65" 

67 

“6-8“ 

69 


Al=(X<b/XNU>**IUIHYD  — 

T  E  M  P  =  1 . 5*C<2*(3.*C0S(XINCL)**2— 1.)/ (1.-E0*E0)**1.5- 
TTETT  T=TE  M  P7  (AT  *  ATT 

A0=A1*(1.-DELl*(.5*T0THRD+DELl*(1.+134./81.*DELl))) 


“D  ETTCT=  TTr«TP7TX 0  *“A'Cn 
XN0DP*XN0/(1.+DEL0) 

I'FTTT  JU'P~I  /~X  NOTXMNPD  A  )  75 E  .  .15~50J  IDbt^l 


“B"STXR“=B“STXR“’dTO  .  *  *  1  B  E XPT7XE" 
TSINCE«TS 


IFITAusI 
IF(I0EEP  .EQ 


1  .AND.  (IEPT  .EQ.  1  .OR.  IEPT 

"5TT0 - 


EQ.  2 


TXT 

71 

rr 

73 

-TV- 

75 


- 1 - .UR.  IhPl  .  tU.  4)  )  b'J  iu 

9  IF(IDEEP  .EQ.  0  .AND.  (IEPT  .EQ.  3  .OR.  IEPT  .EQ.  5)) 
- 1 - GU  TO  880 

10  GO  TO  (21,22,23,24,25),  IEPT 
21  CALL-  SGP(lHAG,ISlNCb) 

GO  TO  60 

2  2  CALL  SGP4 (I  FLAG, I SINlb ) 

GO  TO  60 

23  CALL  SDP4(I~FLAG,  I  SINCE-) 

GO  TO  60 

24  CALL  SbPd(IFLAb,ISlNCE? 

GO  TO  60 

25  C  ALL'  bOP3d  FLAG,  I  S1NCE~) 

60  X  =X  *  X  KMP  E  R / A  E 
- Y  =  Y  *XKHPER/"A’E 

Z=Z*XKMPER/AE 

- X  D0TgXD0T*XKHPER/AE*r|V|NPDA/8640CT^ 

YDOT=YDOT*XKMPER/AE*XMNP DA/36 400. 

Z  DOT=ZDOT *'TK MPER/AE*XMNPDA/~8640 OV 


WRITE (6, 705)  T SINCE, X,Y,Z,XDOT, YD OT,ZDOT 
TSINCEaTSINCE+OELT 

I F(AdS(TSINCE)  .GT.  A8S(TF))  GO  TO  2 


- GO  TO  TO - 

700  FORMATd  1,3F10.0) 


- 9% - 

95 

- 7trt— 

1 

F0RMAT(29x,Dl4.a,1X,3F8.4,/,6'X,t-3.7,F8.4, 

4X,F8.7,I2) 

1  X  ,  2  F  1  1.9,  1  X,F5.5;i'2, 

- 9* - 

97 

- 70?  - 

1 

F0RMAT(18X,D14.8,1X,  FT  0t8T"2  (1X,F6.5,I2),/ 

F7. 7, 2 (IX, F8. 4), IX, FI  1.8) 

*7x*c(1X#F8#4)/iX* 

— rs — 

99 

7tn~ 

704 

F  OR  WA  ~C7vX,  AT; 

FORMAT(1H1,A80,/,1X,A80,//,1X,A4,7H  TSINCE, 

— mo - 

101 

- r 

i 

1  4  A  ,  l(1A,  1  0  A  ,  IMT,  1  0  A  ,  1  Mi,  1  **  A  , 

4HXD0T,13X,4HYD0T,1 3X,4HZD0T,/ /) 

— r&2 - 

103 

7t)5 

706 

F0RMATC7F1  7.8) 

FORMAT( A30) 

—rov — 

105 

707"  "F  OTTMATC  7  9  X  ,  AT) 

930  FORMAT ("SHOULD  USE  DEEP  SPACE  EPHEMERIS") 

— rtrc - 

107 

- 7VCT 

950 

rURMAI  (  SHUULU  USt  NEAR  tAKIH  t“ntntK 1 5  J 

FORMAT ("EPHEMERIS  NU MB E R " , I 2 , "  NOT  LEGAL, 

WILL  SKIP  THIS  CASE") 

— 1“0~8 - 

40  9 

- 

“tfR^TE“66 , 9300 

GO  TO  9 

rro - 

1 1 1 

-  8T0~ 

WRIT  E(6, 940") 

GO  TO  10 

1  T2 - 

- 9OTT 

UR1  TE  (6,950?  I  EP~1 

o 

4 

82 

1 


1 1  3 

GO  TO  2 

- 

11  4 

- ETTD 

• 

83 


4 


1 

- 7 

3 

— *-*■ 

5 

~7T 

7 

7 

9 

— nr 

1 1 

TT 

1  3 

“ TT 

1  5 

TT 

1  7 

— VS 

1  9 
- ZTT 

21 


FUNCTION  ACTAN(SINX^COSX) 

— C0MMON/C2/&E2ftA/PI,Pl02/TMOPI,x3Pl0T 

ACT  AN*0  •' 

- IF  (COSX  .EQ.U'.  )  GO  TO  5 

IF  CCOSX.GT.O.  )  GO  TO  1 

- A  C  T  A  N  =  PT” 

GO  TO  7 

i  it  (s'lNx'rroTir:  )  go  to  3 

I F  (SINX.GT.O.  )  GO  TO  7 

ACTAN=TWOPI 

GO  TO  7 

5  if  (sinx.'E'qTTt;  )  go  to  a 

IF  (SINX.GT.O.  )  GO  TO  6 

- ACT  AN  =  X  3  PI  (52 

GO  TO  8 

6  ACT  AN*P I  02  — 

GO  TO  8 

7  TEMPOS  I NX/COTX 

ACTAN  =  ACTAN+AT  AN(TEMP) 

3  R'ET'URN — — 

END 


FUNCTION  FMOD2PCX) 

COrtHON/C2/DE2ftA,PI,P102,TwOPI/X <Pl02 

F  MO  D  2  P  =  X- 

I =FmOO2P/Tw0PI 

FMOD2P=FMOD2P-I*ThJOPI 
I  F  C  FMOD^P  .  LT  .0.5  FM0D2P=FMQD2P+Td0PI 
RETURN 


1  FUNCTION  THETAG(EP) 

- 2 - COMMON  /  El  /XMO/XNODEO/OMEGAO/EOVx  INcL,xNO,XND  T20,XNDD60,bs  T  AR, 

3  1  X,Y,Z,XDOT,YDOT,ZDOT,EPOCH,D$50 

— — ^  DOUBLE  PRECISION  EPO  CH  *  D  *  T H  E  T  A ,  T  W  OP  I  /  Y R  ,  T  EM  P ,  E  P  / D  S5  0 

5  TW0PI=6. 283185 30 717959D0 

6  YR=CEP+2.D-7)*1.D-3  “ 

'  7  J  Y= YR  _ 

3  Y  R  =  J  Y 

9  D»EP-YR*1.D3  _ _ 

fO  IFCJY.LT. 10)  J  Y= J Y  +  80 

11  N=CJY-69)/4 

- T  2  IFCJY.LT.70)  N  =  CJY-72)A 

13  D  S50a73Q5 . DO  +  365 .DO* C J Y-70 )  +N  +  D 

- r 5 - TH6TA  =  1  .  729444940'0  +  6Y3W3S  809  3 7  DU  *  D  S 50 

15  TEMP=THETA/TWOPI 

- V6 - 3T=TTMP  ' 

17  T  EM Ps I 

T8 - THETAG=rHETA-TEMP*TwOP"r  7 

19  I FCTHETAG.LT. 0. DO)  T HE T AG  =  TH E T A G  + TW 0 P I 

- 2H - RETURN 

21  END 


12.  USERS  GUIDE,  CONSTANTS  AND  SYMBOLS 

The  first  input  card  is  the  program  card.  The  format 
is  as  follows: 


Column 


2-11 

12-21 

22-31 


Format 

II 


F  10.0 
F  10.0 
F  10.0 


Description 

Ephemeris  program  desired 

1  =  SGP 

2  =  SGP4 

3  =  SDP4 

4  =  SGP8 

5  =  SDP8 

Prediction  start  time 
Prediction  stop  time 
Time  increment 


All  times  are  in  minutes  since  epoch  and  can  be  positive  or 
negative.  The  second  and  third  input  cards  consist  of  either 
a  2-card  transmission  or  2-card  G  type  element  set.  Either 
type  can  be  used  with  the  only  condition  being  that  the  two 
cards  must  be  in  the  correct  order.  For  reference  a  format 
sheet  for  the  T-card  and  G-card  element  sets  follows  this 
section. 

The  values  of  the  physical  and  mathematical  constants 
used  in  the  program  are  given  below. 


Variable  name 


CK2 


CK4 


E6A 


Definition 

—  J  a  ^ 

2  2  aE 

-  1  J.  aF4 
4  E 


10 


8 

-6 


Value 


5 . 413080E-4 


62098875E-6 


1.0  E-6 
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Q0MS2T 

(q0  -  s)4  (er)4 

1.8802791 6E -9 

S 

s  (er) 

1.01222928 

TOTHRD 

2/3 

.66666667 

XJ3 

J3 

- .  253881E-5 

XKE 

v  c  er  ->3/2 
ke  CiIH° 

.  743669161E-1 

XKMPER 

kilometers/Earth  radii 

6378.135 

XMNPDA 

time  units/day 

1440.0 

AE 

distance  units/ 

Earth  radii 

1.0 

DE2RA 

radians /degree 

.  174532925E-1 

PI 

IT 

3.14159265 

PI02 

tt/2 

1.57079633 

TWOPI 

2tt 

6.2831853 

X3PI02 

3  tt/2 

4.71238898 

where  er  =  Earth  radii.  Except  for  the  deep-space  models, 
all  ephemeris  models  are  independent  of  units.  Thus, 
units  input  or  output  as  well  as  physical  constants  can 
be  changed  by  making  the  appropriate  changes  in  only  the 
DRIVER  program. 

Following  is  a  list  of  symbols  commonly  used  in  this 
report. 
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n  =  the  SGP  type  "mean"  mean  motion  at  epoch, 
o 

e  =  the  "mean"  eccentricity  at  epoch  v 

o 

i  =  the  "mean"  inclination  at  epoch 
o 

M  =  the  "mean"  mean  anomaly  at  epoch 
o 

w  =  the  "mean"  argument  of  perigee  at  epoch 
o 

n  =  the  "mean"  longitude  of  ascending  node  at  epoch 
o 

n  =  the  time  rate  of  change  of  "mean"  mean 
o 

motion  at  epoch 

n*  =  the  second  time  rate  of  change  of  "mean" 
o 

mean  motion  at  epoch 
B*  =  the  SGP4  type  drag  coefficient 

k  =  Vgm  where  G  is  Newton’s  universal 
e 

gravitational  constant  and  M  is  the  mass 
of  the  Earth 

aP  =  the  equational  radius  of  the  Earth 

.Ci 

J2  =  the  second  gravitational  zonal  harmonic 
of  the  Earth 

J3  =  the  third  gravitational  zonal  harmonic 
of  the  Earth 

J4  =  the  fourth  gravitational  zonal  harmonic 
of  the  Earth 

(t  -  t  )  *  time  since  epoch 
k2  =  7  J2  aE 
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.  3  T  a  4 

k4  *  '  8  J4  aE 


A3,0  J3  aE 

q  *  parameter  for  the  SGP4/SGP8  density  function 
s  *  parameter  for  the  SGP4/SGP8  density  function 

B  »  i  Cn  the  ballistic  coefficient  for  SGP8 
Z  Dm 

where  is  a  dimensionless  drag  coefficient 
and  A  is  the  average  cross-sectional  area  of 
the  satellite  of  mass  m 
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SDC/SCC  TWO/TMREE  CARD  ELEMENT  SET  -  INTERNAL  FORMAT 


Q  1 


I  • 


4 


ELEMENT  SET  FORMATS 


•  4 


13.  SAMPLE  TEST  CASES 


For  reference  a  sample  test  case  is  given  for  each  of 
the  five  models.*  The  input  used  was  standard  T-cards  and 
the  output  is  given  at  360  minute  intervals  in  units  of 
kilometers  and  seconds. 

When  implemented  on  a  given  computer,  the  accuracies 

with  which  the  test  cases  are  duplicated  will  be  dominated 

by  the  accuracy  of  the  epoch  mean  motion.  If,  after  reading 

*  k 

and  converting,  the  epoch  mean  motion  has  an  error  An=jxl0 
radians/t ime ,  then  the  predicted  positions  at  time  t  may 
differ  from  the  test  cases  by  numbers  on  the  order  of 

Ar  =  An (t  -  tQ)  (6,378.135)  kilometers 


*The  test  cases  were  generated  on  a  machine  with  8  digits 
of  accuracy.  After  a  one  day  prediction,  the  test  cases 
have  only  5  to  6  digits  of  accuracy. 
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1  88888U 

2  88888  72.8435 

SGP  TSINCF 

0. 

360. 00000000 
720.00000000 
1080.00000000 
1440.00000000 


80275.98708465  .00073094  13844-3 

1  15.9689  00Q6731  52.6988  l'l0.57l4 


66816-4  0  8 

1  6 .0582-4  5  1  8  105 


2328.96594238 

2456.00610352 

2567.39477539 

2663.03179932 

2742.85470581 


-5995.21600342 

-6071.94232177 

-6112.49725342 

-6115.37414551 

-6079.13580322 


1719.97894287, 

1222.95977784 

713,97710419 

195.73919105 

-328.86091614 


XDOT 

2.91110113 

2.67852119 

2.43952477 

2.19531813 

1.94707947 


YDOT 

-0.98164053 

-0.44705850 

0.09884824 

0.65333930 

1.21346101 


ZDOT 

-7.09049922 

-7.22800565 

-7.31889641 

-7.36169147 

-7.35499924 
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1  8  8  8  8  8  80275.98708465 

2  88888  72.8435  115.9689  0086731 


.00073094  13844-3 

52.6988  110. 571* 


S6P4  TSlNCt 


X 


Y 


* 


0. 

360.00000000 

720,00000000 

1080.00000000 

1440.00000000 


2328.97048951 
2456.10705566 
2567.561g5068 
2663 .09078980 
2742.55133057 


-5995.22076416 

-6071.93853760 

-6112.50384522 

-6115.48229980 

-6079.67144775 


XDOT 

2.91207230 

2.67938992 

2.44024599 

2.19611958 

1.94850229 


YCOT 

-0.98341546 

-0.44829041 

O.O98IO869 

0.65241995 

1.21106251 


66816-4  0  8 

6;058245lF  105 

Z 

1719.97067261 

1222.89727733 

713.96397400 

196.39640427 

-326.38°95856 


ZDOT 

-7.09081703- 

-7.22879231 

-7.31995916 

-7.36282432 

-7.35619372 


1  1 1  an  i  u 

2  1 1  BO  1  46  #  79 1 6 

SDP4  tsince 

0. 

360.00000000 

720.00000000 

1080.00000000 

1440.00000000 


80230.29629788 
230.4354  7310036 

X 

7473.37066650 

-3305.22537232 

14271.29759766 

-9990.05883789 

9787.86975097 


XDOT 

5.10715413 

-1.30113538 

-0.32050445 

-1.01667246 

-1.09425066 


.01431103  00000- 

47.4722  10.4117 

Y 

420.95261765 

32410.86328175 

24110.46411133 

22717*35522461 

33753.34667969 


YOOT 

6.444o8784 

-1.15131518 

2.67964074 

-2.29026759 

0.92358845 


14311-1 

2.2853784b 

Z 

5828.747863  11 
-24697. l76757rtl 
-4725.76837158 
-236l6*t9062501 
-15030.81176753 


ZDOT 

-0.18613096 
-0.28333528 
-2 .08405269 
0.72892364 
-1.52230928 
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4 


1  8  88  88'-'  80275.9870fi465  .00073094  13844- 

2  88888  73.8435  115*9689  0086731  52*6988  i‘10.5714 


SGP8  TSINCE 

0. 

360.00000000 

720.00000000 

1080.00000000 

144Q.OOOOOOOO 


X 

2328.87265015 

2456.04577637 

2567.63383789 

2663.49508667 

2743.29238892 


Y 

-5995.21289063 

-6071.90490772 

-6112.40881348 

-6115.18182373 

-6078.90?83691 


XDOT 

2.91210661 

2.67936245 

2.43992555 

2.19525236 

1*94680957 


YDOT 

-0.98353850 
-0. 44820fl47 
0.09893919 
0.65453661 
1*21500109 


66816-4  0  8 

I6.C58245I8  IO5 

Z 

1720.04684338 

1222.84086609 

713.292a2379 

194.62816810 

-329.73434067 


ZDOT 

-7.09081554 

-7.22888553 

-7.32018769 

-7.36308974 

-7.35625595 


» 
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1  llflOlU 

2  1 l 8° 1  ^ 6 • 79 1 6 

B0230.296297H8 
230,435^  73l0°36 

.01431103  00000-0 

47.4722  10.41 17 

14311-1 

2*2h537848 

SDPO  TSlMCp 

X 

Y 

Z 

o  o  o  o 

O  O  CO 

o  o  o  o 

C  CO  o 

o  o  o  o 
o  o  o  o 
e  o  co 
o  o  o  o 

•  •  •  •  • 

e  c  o  o  o 

ON  cc  vf 
rr,  o  <r 

7469,47631836 
-3337.30992310 
U226. 54333496 
-10151, 59838C67 
9420.08203125 

4 1 5 . 99  3  g07  q2 
32351.39086914 
24236.08740234 

22 223.69848633 
33847.21875000 

5829.64318848 

-24658,63037109 

-4856.19744673 

-23392.39770508 

-15391.06469727 

XDOT 

YDOT 

7  DOT 

5,1  1402285 
-1 .30200730 
-0.33951668 
-1.00112480 
-1.1  1986055 

6.44403201 
-1.15603013 
2.65315416 
-2.33532837 
0.854101 4^ 

-0.18296110 

-0.28164955 

_2.08ll4l53 

0.76987664 

-1.49506933 

# 
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14.  SAMPLE  IMPLEMENTATION 


These  FORTRAN  IV  routines  have  been  implemented  on  a 
Honeywell-6000  series  computer.  This  machine  has  a  pro'- 
cessing  speed  in  the  1MIPS  range  and  a  36  bit  floating 
point  word  providing  8  significant  figures  of  accuracy  in 
single  precision.  The  information  in  the  following  table 
is  provided  to  allow  a  comparison  of  the  relative  size 
and  speed  of  the  different  models*. 

core  used  CPU  time  per  call  (milliseconds) 


Model 

(words) 

Initialize 

Continue 

SGP 

541 

.8 

2.7 

SGP4 

1,041 

1.9 

2.5 

SDP4 

3,095 

5.1 

3.6 

SGP8 

1,601 

00 

• 

ri 

2.2 

SDP8 

3,149 

5.4 

3.2 

« 

r 


*  The  timing  results  are  for  the  test  cases  in  section 
thirteen  with  a  one  day  prediction.  Times  may  vary  slightly 
with  orbital  characteristics  and,  for  deep-space  satellites, 
with  prediction  interval. 
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